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EXECUTIVE  SUMMARY 

PROGRAM  ELEMENTS  AN!)  OVER  AM  OHJECTIVES 


NAP's  current  project,  for  which  this  Semiannual  Report  is  an  interim  report,  has 
thrvc  major  components  and  sets  of  objectives: 

1.  Integrated  1)01)  Voice  and  Data  Networks  - 

The  project’s  specific  concern  is  to  identify  the  appro- 
priate mix  of  switching  technologies  (e.g.,  circuit,  packet, 
and  hybrid  switching)  that  can  best  meet  DOI)  data  and 
voice  communication  requirements  in  separate  or  in- 
tegrated future  networks.  To  perform  the  analyses 
necessary  for  this  determination,  the  project’s  goals  are  to 
identify  key  issues  and  parameters,  to  develop  cost/pcr- 
formancc  tradeoffs  and  technology  assessments,  and  to 
provide  detailed  recommendations  for  best  meeting  pro- 
jected DOD  requirements  for  voice  and  data  communica- 
tions in  the  198(J’s  and  beyond. 

2.  Gatewuy  Topological  Optimization  for  Interconnecting  Packet 
Switched  Networks  - 
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Growth  of  different  packet  switching  networks  within  the 
DOR  i>  leading  to  a number  of  * |u<M D . 1 1 s concerning  the 
most  effective  means  to  connect  these  networks  \ hen 
necessary.  A major  goal  of  thK  project  is  to  develop 
techniques  for  optimizing  the  number  and  locations  of 
gateways  required  to  interconnect  packet  networks. 
Additionally,  these  techniques  will  be  utilized  to  identify 
fundamental  parameters  influencing  gateway  locations, 
and  as  a concrete  example,  gateway  strategies  for 
connecting  ARPANET  and  AUTODIN  II  will  be  recom- 
mended. 

3.  Ground  Packet  Radio  Technology  - 

ARPA  has  developed  a ground  packet  radio  system  that  is  currently  nnd',:,goi?r  a 
sequence  of  experimental  tests  and  performance  verifications.  Tests  are  constrained; 
because  of  the  limited  number  of  repeaters  which  presently  exist  and  because  of  the 
complexity  of  simulating  all  possible  environments  and  stresses  experimentally.  The 
goals  of  this  project  are: 

To  determine  (via  simulation  and  analysis)  the  performance 
profile  of  packet  radio  networks  as  a function  of  key 
system  parameters  such  as  maximum  capacity,  error 
rates,  and  equipment  limitations. 

To  provide  specific  recommendations  for  the  enhancement 
of  the  experimental  packet  radio  system  in  order  to 
improve  factors  such  as  the  time  required  to  stabilize 
after  element  failures  and  system  capacity  under  noisy 
conditions. 
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To  determine  the  speed  mid  eharaeterist ics  of  device's 
within  a packet  radio  system  whose  elements  are*  mobile. 

To  determine  the  range  for  which  the  pns  at  oetwoik 
procedures  hi  come  infeasible  or  introduce  unacceptably 
high  performance  degradations,  and  to  prcpjse  feasible 
system  design  alternatives  for  use  within  mobile  networks. 

This  document  contains  only  those  interim  results  which  arc  both  complete  and 
of  utility  ns  stand  alone  items.  At  the  completion  of  the  contract,  these  results  will  be 
intern  'ted  into  the  filial  report. 

Results  given  in  this  interim  document  are  reported  in  stand  alone  chapters  and 
arc  related  to  one  or  more  of  the  major  program  elements  with  their  associated  lasks. 
Chapter  1,  ,r Topological  Design  of  Mixed  Media  Networks,”  is  relevant  to  both  the 
Integrated  Voice  and  Data  and  the  (iateway  Projects.  Chapter  2,  "An  AJ  oritlnn  for 
Design  of  Non-llierarchicul  Circuit-Switched  Networks,”  m a major  ingredient  of  tin* 
Integrated  \oiee  and  Data  Project.  Chapters  3 and  4 arv  part  of  the  Packet  Radio 
effort:  Chapter  3.  ”On  Connectivity  in  Mobile  Packet  Radio  Networks,”  relates  to 

system  mobility  issues;  Chapter  4,  ”An  Approximate  Analytical  Model  for  Initialization 
of  Single  Hop  Packet  Radio  Networks,”  is  part  of  the  effort  to  develop  a performance 
profile  for  the  Packet  Radio  System. 

The  project  plan  for  each  of  the  three  program  elements  is  directed  towards: 
identifying  fundamental  issues  and  problems  related  to  overall  objectives;  structuring 
specific  tasks  whose  completion  will  solve  these  problems;  accomplishing  the  tasks; 
and  integrating  the  results  of  the  tasks  to  meet  the  overall  program  objectives.  Major 
tasks  identified  are  given  below: 

1.  Integrated  DOL)  Voice  and  Data  Network  Tasks 

-1.1  Develop  performance  measures  for  packet  voice  and  determine 
the  impact  on  switch  architecture  (hardware  and  software). 

-1.2  Develop  techniques  for  the  design  and  analysis  of  packet 
switched  networks  for  voice  ' pd  integrated  voice  and  data. 
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-1.3  Investigate  the  inipm'l  of  priority  structures  for  voioo  and  data 
in  packet  s v\  i to  I led  networks. 

-1.4  Apply  teehna, lies  developed  nbo\e  In  AlJ'IOWhN  ,uu]  At  TOPIN’ 
II  truffle*  and  generate  performance  profiles  lor  the  imcm-l  switching 
technology. 

-1.5  Develop  approaches  for  the  design  and  analysis  of  circuit 

switching  networks. 

-1.6  Develop  a methodology  for  assessment  of  circuit  switch 

technology  with  regard  to  hardware,  software,  cost  and  capacity. 

-1.7  Apply  the  circuit  switching'  design  methodolo  ;v  to  DO!)  voice, 
data,  and  combined  voice  and  data  require?!. mil  , .md  • in  i ip 
performance  profiles  for  the  circuit  swii  ‘leu  • t«  ehn<  I*  . . 

-1.8  Develop  algorithms  for  the  design  and  analy  T of  integrated 
cireuit/packct  switching  networks.  (Such  mixed  switching  strategies 
are  called  hybrid  strategics.) 

-1.9  Conduct  an  assessment  of  integrated  switches  being  developed, 
in  terns  of  cost,  modularity,  reliability,  and  impact  of  priority 
structures. 

-1.10  Determine  the  cost-effectiveness  of  incorporating  satellite 
subnetworks  into  integrated  packet/circuit  switched  networks. 

-1.11  Determine  partition  criteria  for  classes  of  voice  and  data  to  be 
served  by  the  circuit  and  pucket  switching  components  of  a hybrid 
switching  network. 

-1.12  Apply  algorithms  to  AUTOVON  and  AUTODIN  II  data  bases  and 
generate  performance  profiles  for  the  hybrid  switching  technology. 
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-1.13  Integrate  the  findings  of  Tasks  1.4,  1.7.  mid  1.12  mid  provide 
recommendations,  sensitivity  studies,  and  cost  and  performance 
comparisons  between  tic  candidate  packet,  circuit,  and  hybrid 
alternatives. 


2.  Gateway  Topoiojpival  Optimization  Tasks 

-2.1  Determine  issues,  parameters,  and  performance  criteria  for 
interconnecting  packet-switched  networks. 

-2.2  Develop  a methodology  and  computer  programs  for  determining 
the  number  and  location  of  gateways  for  interconnecting  either 
terrestrial  or  terrestrial  and  satellite  networks. 

-2.3  Apply  the  methodology  to  a ease  sturv  - the  interconnection  of 
ARPANET  and  \UTuDIN  II. 

3.  Ground  Packet  Radio  Technology  Tasks 

-3.1  Determine  the  performance  and  capacity  profile  of  packet  radio 
networks  for  different  error  rates,  code  rales  and  topologies  via 
simulation  techniques. 

-3.2  Develop  analytical  models  and  study  network  initialization  time 
ns  a function  of  system  parameters. 

-3.3  Compare  analytical  results  with  simulation. 

-3.4  Upgrade  the  packet  radio  simulator  to  contain  functions  and 
capabilities  such  as  error  correction,  protocols,  and  packet  handling 
techniques  being  implemented  in  the  experimental  system. 

-3.5  Perform  simulation  studies  for  the  elements  included  in  Task 
3.4  above  parallel  to  experimental  tests  and  provide  guidelines  for 
enhancement  (e.g.,  capacity  increases  or  delay  reduction)  of  the 
experimental  packet  radio  system. 
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-3.fi  Determine  the*  speed  and  mobility  characteristics  of  devieex  in 
tin*  packet  radio  system  fur  wind)  the  present  routing  and  miti.ih/a- 
tion  algorithms  degrade  performance  or  become  infeasible. 

-3.7  Propose  rooting  and  iniliali  ation  procedures  for  mobile  packet 
radio  networks  which  will  mitigate  degradations. 

RESULTS 

CHAPTER  I:  TO  PO  LOCI  1C  A L DESK IN  OF  MIXED  MEDIA  NETWORKS: 
DETERMINATION  OF  SATELLITE  AND  TERRESTRIAL  BACKBONE  NODES 

A computer  methodology  is  described  for  optimizing  the  location  of  [racket 
switches  and  satellite  ground  stations  in  a mixed  terrestrial  satellite  MOD  network. 
The  techniques  described,  which  have  been  programmed  and  lorded  for  the  examples 
discussed  below,  combine  heuristic,  comimiatori.il.  and  analytic  elements.  11, e 

algorithms  developed  have  direct  relevance  to  ti  e general  gateway  location  problem 
and  will  be  used  in  forthcoming  efforts  for  this  problem.  Results  reported  here 
include: 


Basic  parameters  which  influence  cost  and  performance  in  a 
mixed  network  are  identified.  Critical  parameters  include: 
Ratios  of  shortest  path  lengths  within  network  to  direct 
distances,  line  overhead,  average  link  utilization,  unit  terres- 
trial and  satellite  channel  costs,  and  average  channel  length.  A 
fundamental  relationship  between  cost  and  traffic  requirements 
is  developed  which  can  provide  network  cost  estimates,  obtain- 
ed without  detailed  layout,  to  within  10%  of  the  best  design. 
This  relationship,  which  is  an  analytic  combination  of  the  basic 
parameters,  is  the  first  such  result  for  packet-switched  sys- 
tems. 


The  techniques  developed  are  applied  to  the  AUTODIN  II  data 
network  problem.  Results  include: 
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Network  cost  (excluding  operations,  encryption  mui  other 
add  on  cost  lectors)  without  satellites  is  on  the  order  of 
$(.*.2  million  per  year.  This  is  a significant  fraction  of  the 
AUTODIN  II  cost. 

Satellites  used  in  backbone  network  reduce  communica- 
tions cost  by  about  7%  (i.e.,  over  $0.6  million  per  year)  for 
DCA  anticipated  user  traffic  requirements.  This  relation- 
ship is  expected  to  be  true  for  other  similar  applications. 

Critical  cost  elements  which  influence  cost  tradeoffs  are 
terrestrial  line  and  satellite  ground  station  costs.  The 
satellite  space  segment  and  other  costs  do  n >t  have 
significant  impact.  Therefore,  further  reductions  over  7 v, 
for  the  given  data  levels  re  achievable  primarily  thi.-u  ,h 
reemeod  cost  ground  stations  and/or  lower  tariffs. 

For  this  application,  the  best  number  of  satellite  ground 
stations  is  sm  11  and  constant  (at  4)  as  the  number  of 
backbone  picket  switches  increase  from  7 to  16.  We 
expect  that  this  will  be  true  for  other  applications  unless 
ground  station  cost  is  significantly  lower  (e.g.,  1/5  to  1/10 
of  current  costs). 

For  this  application,  and  for  designs  with  the  same  number 
of  backbone  nodes,  those  with  lower  local  access  cost 
nearly  always  have  lower  total  cost  (i.e.,  minimum  local 
access  cost  nearly  always  implies  minimum  total  network 
cost). 
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CHAPTER  2:  AN  M.dOKITIIf!  FOR  PI-SICK  OF  NONIlIFU  ARCIIK ' \l, 

circui  r .switched  n ktworks 


\ general  eomput  ri/otl  procedure  tor  the  minimum  »*ost  design  of  circuit 
switched  networks  is  reported  and  studies  to  he  conducted  with  the  procedure  me 
described.  I he  procedure,  which  was  developed  to  enable  the  performance  of  these 
studies,  needed  because  no  organization  has  previously  required  efficient  tools  witn 
which  to  examine  all  of  the  basic  assumptions  involved  in  trndilioial  circuit  switched 
network  design  unfctcred  by  the  constraints  that  have  evolved  piecemeal  for  existing 
circuit  switched  systems  such  as  the  U.  S.  telephone  network.  In  this  chapter  we 
report  on  the  design  of  the  procedure.  Future  studies  with  the  system  will  identify 
cost  performance  tradeoffs  for  voice  and  data  as  a function  of  traffic  load  and  mix, 
routing  strategy,  and  structure  of  the  circuit  switch  (switch  set  up  time,  signalling 
scheme,  etc.) 

Relevant  features  of  the  procedure  are: 

Computationally  elhcienl  and  thus  useful  for  a wide  range  of 
studies,  including  investigation  of  switch  design,  signalling 
scheme,  and  network  layout. 

Inputs  are  switch  locations,  traffic  requirements,  and  system 
constraints  such  as  available  tariffs  and  desired  user  call 
rejection  probability. 

Outputs  are  network  link  layout,  link  capacity  assignment, 
routing  plan  and  network  cost  to  meet  traffic  requirements  and 
constraints. 

Contains  separate  modules  to  enable  performance  analysis, 
alternate  routing,  trunk  sizing,  and  line  layout.  Thus,  useful  for 
studying  impact  of  variations  in  switch  and  network  operating 
procedures. 
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CilAPTFK  3:  ON  ( 'ON  N \'A  ' 1 IV1TY  IN  MOHILK  PACKKT  KADIO  NI-TTVOKKS 


This  chapter  reports  the  first  results  on  Packet  Kadio  System  performance  v. it.; 
mobile  (Menu  nts.  An  important  issue  for  sueli  networks  is  lie*  : 1 1 1 • ' > 1 n 1 1 of  ovorho.-wi  data 
that  must  be  sent  through  the  network  for  the  routing  and  eontrol  procedures  to 
operate.  Speed  beeoines  a faetor  that  can  degrade  network  performance  because  of 
this  overhead.  Performance  can  eventually  deteriorate  to  the  point  where  operations 
become  infeasible.  Consequently,  it  is  important  to  design  networks  which  exhibit 
graceful  degradation  over  as  wide  a range  ns  possible. 

The  techniques  described  in  this  chapter  will  be  used  in  forthcoming  efforts  to 
study  performance  of  current  and  proposed  network  operating  rules  to  identify 
mobility  factors  leading  to  infeasibili tv.  Procedures  are  developed  to  calculate  the 
following  quantities  in  systems,  whose  element  trajectories  are  known  or  can  be 
computed: 


rimes  during  which  a repeater  can  communicate  with  the 
station. 

Time  interval  over  which  the  entire  network  is  connected  (i.e., 
all  repeaters  have  paths  to  the  station). 

Connectivity  profile  of  the  network  (time  spans  when  "at  most," 

"at  least"  or  "exactly"  a specified  number  of  repeaters  are 
unable  to  communicate  with  station). 

The  trajectory  computability  assumption  will  he  relaxed  in  future*  studies  and  the 
performance  projections  currently  being  calculated  will  serve  on  bounds  on  system 
performance. 
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CJ1IAT  i ICR  -I:  AN  APPROMM  AT  I-  AN  AFYTICA  I.  .V,Ol)i:L  I OK  INI  IT  A UZATTON 

ofsinou:  nor  paokkt  radio  ni  t works 

This  chapter  deals  with  initiaiiza lion  in  Packet  Radio  Systems  mu!  is  part  of  the 
effort  to  develop  a performance  profile  of  the  system.  An  analytic  model  is  first 
developed  to  calculate  initialization  times  (the  time  required  to  stabilize  the  network 
given  a disturbance  such  as  the  entry  or  failure  of  repeaters).  The  model  is  then  used 
io  identify  efficient  operating  pnrnmctc.s  such  as  those  listeo  below  for  the  Packet 
Radio  System.  Future  effopts  will  extend  the  model  to  general  packet  radio  systems 
and  be  used  to  investigate  the  effect  of  initialization  procedures  on  the  time  required 
for  network  stabilization  under  various  operating  conditions.  Current  numerical 
results  include  operating  parameters  for: 

Station  transmissions  of  connectivity  information  (rv;..  a libel 
transmission  every  forr  maximum  packet  times  i mo:  l <■  ffee- 
tive). 


Repeater  initialization  packets  (called  KOP’s)  (o.g.,  if  tne 
network  has  M repeaters,  a repeater  slioi  Id  transmit  a HOP 
every  m(e)  maximum  packet  times). 
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CHAPTER  1 

topoi, on i cal  deb i gn  op  mixed  media  networks : ni:rr i;rm i nation 

OF  SATELLITE  AND  TERRESTRIAL  BACKBONE  NODES 


1 . INTRODUCTION 


In  this  chapter,  we  address  the  topological  design  of  large, 
distributed  data  net vorks  m which  both  terrestrial  and  satellite 
technologies  arc  employed  - which  we  called  mixed  media  data  net- 
works. The  particular  problem  considered  is  the.  determination  of 
terrestrial  and  satellite  backbone  switch  locations  which  minimise 
total  network  cost. 

In  the  following,  wo  first  briefly  consider  the  various 
options  in  packet  switched  satellite  communication.  We  then 
consider  the  general  topological  design  problem  of  large  distrib- 
uted networks,  and  indicate  how  the  switch  location  problem  ties 
into  the  general  design  problem. 

1 . 1 Packet-Switched  Sate Hite  Communi cation 

Until  now,  the  packet-switched  networks  implemented  use  mainly 
terrestrial  communications  links.  However,  with  the  advancement  of 
satellite  communication  technology,  -the  possibility  of  using  satel- 
lite links  to  reduce  total  network  cost:  and  increase  network  growth 
flexibility  has  been  receiving  increasing  attention  [ABRAMSON,  1973] , 
[GERLA,  1974b),  [IIUYNII,  1976],  [KLEINROCK,  1.973],  [LAM,  1974], 
[ROBERTS,  1973). 

Satellite  communications  can  be  used  in  the  following  modes  in 
a packet-switched  environment.  (Here,  we  only  consider  alternatives 
at  the  packet-switched  backbone  level.) 
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1.  Dedicated  poi  i t- 1 .o-po  ; n L connect  ionr.  between  earth 
stations  at  some  of  the  backbone  node  sites,  rcpl.icincj 
more  costly  terrestrial  clmnneln.  This  mock.'  is  not  very 
attractive  unless  there  are  very  few  earth  stations  (say 
2 or  3 ) . 

2.  Dedicated  up-links  from  earth  stations  located  at 
some  of  the  backbone  node  sites  to  the  satellite,  and 
broadcast  down-links  lrom  the  satellite  to  all  earth  sta- 
tions . 

3.  Multiple-access  (ALOHA # slotted  ALOHA,  reservation) 
up-link  from  the  earth  stations  to  the  satellite  channel, 
and  broadcast  down-links  from  the  satellite  to  the  earth 
stations . This  mode  is  attractive  when  the  users  (earth 
stations)  have  different  busy  hour,  or  when  there  is  a 
cost-effective  high  bandwidth  satellite  channel  option. 

1 . 2 Design  Approaches  for  barge  Distributed  Data  Networks 

Large  distributed  networks  generally  result  from  the  integra- 
tion of  existing  centralized  and/or  distributed  data  communications 
systems.  The  purpose  of  integration  is  to  achieve  line  economies, 
intersystem  communications  capability,  higher  network  bandwidth, 
improved  reliability,  improved  growth  flexibility,  and  more  effi- 
cient network  control  and  management  [MAC,  1976].  Consolidation 
and  integration  of  communication  requirements  is  presently  being 
considered  by  several  large  corporations  as  well  as  government  and 
military  agencies.  For  example,  the  most  substantial  packet  net- 
work under  development  is  the  AUTODIN  II  System  whereby  the  Depart- 
ment of  Defense  is  in  the  process  of  integrating  many  of  its  ADP 
systems  into  a large,  hierarchical  network. 
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Although  specific  network  i inpl  omental  i on.f;  vary  conn  i durably 
do  pond  i iuj  on  the  application,  Lho  conunon  structure  of  large  distri- 
buted networks  is  the  multi  level  hierarchical  structure  with  a back 
bone  network  at  the  higher  level,  and  local  access  networks  at  the 
lower  levels.  Backbone  and  local  access  nets  may  be  further  sub- 
divided into  hierarchical  sublevels. 

The  backbone  network  is  characterised  by  distributed  traffic 
requirements  and  is  generally  implemented  using  the  packet-switch- 
ing technology.  Local  access  networks,  on  the  other  hand,  have  a 
centralized  traffic  pattern  (all  the  traffic  is  to  and  from  the 
gateway  backbone  node)  and  are,  therefore,  implemented  with  conven- 
tional teleprocessing  techniques  such  as  multiplexing,  concentra- 
tion, and  polling. 

The  selection  of  the  most  effective  network  architecture  (i.e. 
number  of  levels  and  type  of  access  at  each  level)  is  the  first 
problem  that  must  be  attacked  in  the  design  of  a large  network.  A 
discussion  on  alternative  architectures  can  be  found  in  [NAC,  1976] 

Having  selected  the  architecture,  we  must  then  design  a cost- 
effective  topology  that  minimizes  lire  and  nodal  processor  costs 
within  such  architecture.  For  the  multilevel  case,  the  problem 
may  be  subdivided  into  four  subproblems: 

a.  Preliminary  clustering  of  the  user  installations; 

b.  Selection  of  backbone  nodal  processor  locations; 

c.  Local  access  design; 

d.  Backbone  topoloyy  design  for  the  higher  level  net- 
work . 
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One  approach  for  the  i opol  oy  ica  1 cicr.  ign  of  mixed  media  network 
is,  in  subproblem  b of  the  above,  determining  not  only  the  backbone 
nodal  processor  locations,  but  at  the  same  Lime  which  of  them  are 
satellite  switch  locations.  However,  as  indicated  in  Section  3.2, 
the  tradeoffs  involved  here  are  quite  complex.  Thus,  we  are  con- 
tent with  only  determining  the  backbone  switch  locations  in  subprob 
lem  b,  and  then  selecting  the  satellite  switches  from  the  sot  of 
backbone  switches  in  subprobl  cm  d.  Consequently,  subproblem  d is 
further  broken  into  the  following  two  parts: 

d.l  Selection  of  satellite  switch  locations; 

d.2  Topology  design  of  the  satellite  subnet 
(satellite  channel  and  earth  stations)  and 
terrestrial  subnet. 


Subproblem  a,  the  partitioning  of  a large  terminal  population 
into  "minimal  cost"  clusters  satisfying  given  constraints  is  a 
classical  problem  in  data  network  design.  An  efficient  technique 
for  clustering  a population  of  multidropped  terminals  was  present- 
ed in  [MCGREGOR,  197  5]  . Similar  techniques  may  be  constructed 
for  other  local  access  strategies. 

Subproblcm  c,  the  design  of  local  access  topologies,  has  also 
been  well-studied.  Several  algorithms  are  found  in  the  litera- 
ture [CHOU,  1973],  [ESAU,  I960],  [WOO,  1973].  Typically,  the 
algorithms  solve  the  problem  of  optimal  location  of  concentrators, 
and  optimal  layout  of  multidrop  lines  to  connect;  terminals  to 
concentrators . 

Subproblem  c,  the  selection  of  the  number  and  location  of 
backbone  switches  has  recently  been  studied  extensively  lor  the 
terrestrial  network  design  in  [NAC,  1976].  Two  algorithms:  ADD 

and  CLUSTER,  are  proposed.  xn  this  chapter,  we  propose  a simpli- 
fied backbone  cost  estimate  for  the  ADD  algorithm  (Section  3.2), 
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and  apply  l ho  simplified  ADD  algorithm  to  the  baekbonc  swi ‘ h 
selection  problem  for  mixed  media  networks.  It  is  experimentally 
demonstrated  that  the  cost  function  proposed  is  a fairly  good 
approximation  to  that  used  in  [MAC,  J976],  and  the  saving  in 
computation  time  is  significant. 

Efficient  heuristics  for  the  topological  design  of  terrestrial 
backbone  networks  have  been  reported  in  [FRANK,  1972],  [GERLA,  1974a], 
[hAVIA,  1975].  Once  the  satellite  switch  locations  have  been  deter- 
mined, these  procedures  can  be  adopted  to  subproblem  d.2,  the  de- 
sign of  mixed  media  backbone  networks  [GERLA,  1974b]. 

In  this  chapter,  a simple  cost  function  which  formulates  sub- 
problem d.l  and  an  efficient  exhaustive  search  algorithm  are  devel- 
oped. Experimental  results  demonstrate  that  the  cost  function  used 
.is  a good  approximation  to  the  total  network  cost.  Extensive  ex- 
perimental results  which  study  total  network  cost  as  a function  of 
number  of  earth  stations,  traffic  levels,  and  sensitivity  to  cost 
variations  are  also  presented. 

The  work  presented  in  this  chapter  is  within  the  effort  of  the 
general  topological  design  of  interconnecting  communication  networks. 
One  of  the  subproblems  in  the  latter  is  to  select,  in  each  network 
to  be  interconnected,  a subset  of  nodes  for  "gateway"  nodes.  It  is 
expected  that  the  algorithm  presented  for  solving  subproblems  b 
and  d.l  will  be  used  for  the  gateway  selection  problem,  using  a 
different  cost  function. 
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2.  REVIEW  OF  TERRHSTRl AL  BACKBONE  SWITCH  LOCATION  ALGORITHMS 

In  the  previous  phase  of  AREA  contract  [NAC,  1976],  we  have 
investigated  the  switch  location  problem  for  terrestrial  packet- 
switched  networks.  Briefly  stated,  the  problem  is;  given  user 
traffic  requirements,  facility  cost  functions,  and  candidate  sites 
for  backbone  switches,  determine  the  number  and  local. ion  of  back- 
bone switches  so  that  with  appropriate  optimum  backbone  link  topol- 
ogy, the  overall  communication  cost  is  minimized  while  satisfying 
constraints  on  backbone  delay,  switch  capacity,  etc.  The  basic 
underlying  assumptions  are  as  follows: 

1.  The  user  sites  are  proclustorod , so  that  each  site1 
is  connected  to  a backbone  switch  directly. 

2.  Only  one  type  of  backbone  switch  is  used. 

3.  Only  one  type  of  backbone  trunk  is  used  (though 

several  trunks  in  parallel  may  connect  the  same  two 
swi tches ) . 

Two  solution  procedures,  MODULARIZED  ADD  and  CLUATER,  have 
been  proposed  in  (NAC,  1976] . Of  the  two,  ADD  is  found  to  consis- 
tently give  better  results,  though  it  is  also  more  time-consuming. 
The  following  guidelines  apply  in  selecting  the  most  effective  algo- 
rithm: if  the  number  of  candidates  is  much  smaller  than  the  number 

of  user  sites,  the  ADD  algorithm  is  the  best  choice;  if,  on  the 

other  hand,  the  number  of  candidate  sites  is  very  large  (e.g.,  the 
candidate  site  set  consists  of  user  sites) , the  CLUSTER  algorithm 
is  prefer rable. 
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2 1 Backbone'  Net  work  Cost  Estimation 

| Since  the  optimal  node  location  strategy  is  the  result  of  the 

tradeoff  between  local  access  cost  and  backbone  cos L , the  backbone 
topology  should  in  principle  be  redesigned  at  eacli  iteration  to 
evaluate  the  backbone  cost..  This  approach,  however,  is  computation- 
ally too  time-consuming,  especially  if  the  switch  location  algo- 
rithm requires  at  each  iteration  the  design  and  comparison  of  several 
backbone  configurations,  one  for  each  candidate  node  selection.  We 
need,  therefore,  an  approximate  cost  estimate  which  is  both  compu- 
tationally efficient  and  consistent  (in  the  sense  that  the  error 
introduced  is  "systematic",  without  severe  jumps  corresponding  to 
perturbations  in  switch  number  and  location).  The  estimate  need 
not  be  very  accurate  in  t he  absolute  sense,  since  we  are  consider- 
ing only  cost  variations  relative  to  insertion  or  removal  of  one 
no le  at  a time.  This  estimate,  and  the  procedure  used  to  obtain 
it , can  play  a major  role  in  the  overall  problem  of  the  character- 
ization of  networks  via  simple  parameters  for  use  in  the  i n tor net- 
work gateway  1 ocation  problem . 

We  first  introduce  the  concept  of  nondirect  routing  penalty  P. 
The  route  R used  by  the  traffic  from  i to  j depends  on  the  topology, 
the  traffic  conditions,  etc.  However,  we  would  expect  the  mileage 
on  R to  be  proportional  to  the  distance  between  1 and  j,  and,  there- 
fore, may  approximate  the  length  | R | a s : 

i 

* |R|  - P x d..  (1) 

ft  J 

Where  P is  the  "nondirect  routing  penalty"  defined  as  the  ratio 
between  the  shortest  path  distance  and  the  direct  distance,  for  the 
average  source-destination  pair  (P>1) , and  d^  is  the  direct  dis- 
tance from  i to  j . 
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With  this  estimate,  the  backbone  line  cost  D can  be  expressed 

by: 


cr  ■ .d.  . (lib)  I* 

i)  >;  >;  e u + nn  x f 

i i p 

where : 

p = Avera  cio  Jink  utilization  ( p < 1), 


(2) 


F = Fixed  cost  per  node  ( averaqo  number  of  line 

terminations  per  node  x line  termination  cost) , 


c = Cost/mile  x unit  bandwidth  x month, 
r^j  = Traffic  requirement  from  i to  j, 
b = Lino  overhead  (protocols,  etc. ) , [0<b<J ) , 

NN  = Number  of  switches. 


The  coefficients  P,  P and  F are  determined  experimentally  for 
various  values  of  NN.  Application  results  show  that  such  coeffi- 
cients depend  solely  on  NN,  and  are  insensitive  to  switch  reloca- 
tions [NAC,  1976].  Moreover,  experimental  results  [NAC,  1976]  show 
that  the  estimate,  Eq . (2) , is  adequate  for  the  purpose  of  switch  site 
optimization . 

In  the  switch  location  problem  for  the  mixuu -media  networks, 
due  to  the  added  complexity  of  the  backbone  transmission  cost  trade- 
offs, instead  of  the  estimate  d veloped  in  this  Section,  a much 
simpler  estimate  is  adopted  (see  Section  3.2).  However,  the  trans- 
mission cost  estimate  developed  here  is  used  in  the  selection  of 
satellite  switches  (Sectijn  4). 
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2.2  ADD  Approach 

In  the  ADD  approacn,  a " f igura  t ive"  switch  (cen  ler-of-ma  ss  ) is 
created  (and  dynamically  relocated),  such  that  each  user  that  is  not 
yet  assigned  to  any  selected  switches  is  assigned  to  this  node. 

Each  remaining  switch  candidate  location  is  then  evaluated  by  de- 
termining the  cost  reduction  which  would  be  achieved  by  placing  a 
switch  at  the  location  and  profitably  assigning  some  of  the  unas- 
signed users  to  that  switch.  The  cost  tradeoff  is  based  on  Lnc 
throe  cost  components:  backbone  line  cost  (which  is  estimated  by 

Eg.  (2)),  local  access  cost,  and  switch  cost.  The  location  giving 
the  greatest  cost  reduction  is  then  selected  as  the  next  switch 
site,  and  the  user  sites  contributing  to  ics  selection  are  assigned 
to  it.  When  no  further  cost  rcduct.  a can  be  achieved,  the  process 
halts,  and  the  switch  candidate  location  closest  to  trie  "figurative" 
switch  is  selected  as  the  location  for  the  last  switch. 

The  ADD  approach  can  be  improved  (at  tne  expense  of  increased 
computation  time)  by  considering  each  switch  as  composed  of  several 
modules.  During  each  iteration  of  the  ADD  procedure,  instead  of 
adding  one  nioie  switch,  only  one  more  module  is  added  to  the  back- 
bone. We  call  this  the  MODULARIZED  ADD  procedure.  The  module*  cost 
can  be  made  such  that  the  i-th  module  at  a location  costs  only  a 
fraction  of  the  (i-l)-th  module  at  the  same  location.  This  modi- 
fication is  particularly  appropriate  when  the  backbone  switches  are 
very  powerful  but  roL  ively  inexpensive  as  compared  to  the  other 
network  components.  Experimental  results  [NAC,  llJ76]  indicate  that 
the  MODULARIZED  ADD  approach  does  generate  much  better  results. 

Below  we  develop  in  some  detail  the  cost  trade  off  evaluation 
for  selecting  each  new  switch.  Suppose  k candidates  C^#...,C^  have 
already  been  chosjn  as  backbone  switches  (k  may  be  0,  in  which  case 
no  switch  has  been  chosen  yet). 
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! 


II  all  user  locations  are  assigned  to  one  of  the  C^'s,  i ~ 

1, . . . ,k,  then  we  are  done.  Suppose  on  the  other  hand,  that  not  all 
the  user  locations  are  assigned  Let, 

= Set  of  user  locations  not  yet  assigned. 

If  we  select  only  one  more  backbone  switch,  then  a reasonable  choice 
for  this  node,  Cg,  is  the  traffic-weighted  center  of  mass  of  the  nodes 

in  *0- 

Let  Q be  a backbone  candidate  not  yet  used.  If  wo  use  Q as 
Ck-fl'  t*ien  we  would  have  to  take  some  of  the  nodes  in  Cg  (nodes  not 
yet  assigned)  away  from  Cg,  and  assign  them  to  o.  Thus,  rhe  basic 
tradeoff  is  between  homing  nodes  to  Q,  or  homing  nodes  to  Cg. 

Lot  T be  a user  site  in  Cg.  The  saving  oi  homing  T to  Q,  S^(Q), 
is  given  by: 

St(Q)  = LASt(Q)  + BLSt(Q) 
where; 

LASt(Q)  = Local  access  saving  of  homing  T to  Q. 

BLSt(Q)  = Backbone  line  saving  of  homing  T to  Q. 

Tae  local  access  saving  is  given  by: 

LASt(Q)  = Kt  x (d(T,Cg)  - d (T,Q) ] 


where; 

Kt  --  Mileage  cost  of  the  local  access  line  connecting 
T to  the  backbone. 
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The  backbone  line 
other  simplifying 

BLST(Q)  tT 


saving  can  be  estimated,  with  Kgs.  (1)  - 
assumptions  lNAC;  1976],  by: 

CB  k 

XT — - — X >;  t.  X fd  (c  ,c. ) -cl  (c>,c. ) ] 
Total  i=l 


(2)  and 


where ; 


t 


T 


= Total  traffic  (sum  of  transmit  and  receive 
traffic)  at  user  site  T. 


tj  = Total  traffic  of  users  assigned  to  switch 

i 1 , . • • , k • 

t-,  , - Total  traffic. 

Total 


= Cost  factor  as  defined  by  the  approxa ma tion 
method  (C„  = c x P/p,  where  e,  P and  p are  as 

L> 

defined  in  Section  3.2). 


We  proceed  by  assigning  to  Q the  user  sites  with  positive  sav- 
ings, starting  from  the  site  with  the  largest  saving,  until  the 
switch  capacity  is  fully  utilized.  The  saving  of  selecting  Q as 
Ck+1'  t^en  given  by; 

S(Q)  = IS  St(Q)]-C, 

T assigned 
to  Q 

Where  is  the  fixed  cost  of  adding  a backbone  switch. 

After  evaluating  S (Q)  lor  all  the  switch  candidates  0/  the 
candidate  with  the  largest  positive  saving  is  selected  as  the  (k+1) 
switch,  the  process  continues  until  no  candidate  can  give  a positive 
saving  or  if  all  users  are  assigned.  In  this  case,  the  candidate 
nearest  to  the  center  Cr  is  selected  as  the  last  switch  location. 
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2.  3 CMISTF.R  Appronrjh 

The  CLUSTFR  approach  works  Lo  dcLcrmino,  for  a given  va  Luo  of 
N,  the  best  set  of  N locations  at  which  to  place  the  switches.  The 
basic  underlying  assumption  is  that  the  backbone  switch  locations 
that  minimize  local  access  cost  also  minimize  total  network  cost. 
Based  on  the  experimental  evidence  reported  in  Section  5.2,  this 
assumption  is  quite  valid.  The  switch  location  problem  is  thus 
reduced  to  one  of  optimally  partitioning  the  user  sites  into  N sub- 
sets. This  is  hour istical ly  accomplished  by  forming  larger  and 
larger  clusters.  The  clusters  are  formed  by  "rolling  snowballs"  in 
a rather  "balanced"  fashion.  First,  the  two  nodes  closest  together 
are  selected.  These  nodes  are  then  replaced  by  a single  node  at 
their  "center-of -mass"  with  the  combined  weighting  factors  of  the 
first  two  nodes.  The  merging  process  continues  on  a closest  node 
pair  basis  until  only  N nodes  remain. 

in  practice,  the  approach  outlined  above  is  too  simple  to  pro- 
duce a set  of  N good  locations  for  the  backbone  switches.  The  pri- 
mary deficiency  is  the  possible  dramatic  difference  in  size  of  the 
clusters,  which  in  practical  situations  usually  results  in  poor  de- 
signs. To  compensate  for  this  tendency,  a parameter  Z is  defined  as 
a cluster  capacity,  and  a parameter  a is  defined  as  a size  threshold 
expressed  as  a percent  of  Z.  Clusters  will  be  grown  until  they 
reach  a size  greater  than  aZ  and  then  they  will  be  stopped.  Clusters 
are  not  permitted  to  merge  if  the  merger  exceeds  Z.  The  first  N 
cluster  to  be  stopped  will  be  selected  as  the  backbone  nodes. 

By  not  permitting  the  size  Z to  be  exceeded,  the  process  of  de- 
veloping clusters  may  stop  prematurely  for  lack  of  feasible  mergers. 
In  this  case,  simply  the  N largest  clusters  are  used.  Furthermore, 
mergers  may  occur  over  extraordinary  distances  because  of  feasibility 
issues.  This  is  certainly  not  desirable,  and  so  a parameter  of 
maximum  distance  between  mergeable  nodes  is  introduced.  Finally,  if 
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modular  capacity  switcher,  arc  available,  then  several  clufitors  in 
the  immediate  vacinity  of  one  another  art*  probably  hotter  served  by 
one  lnrcje  switch.  Consequently,  a parameter  is  available  to  define 
a minimal  separation  between  clusters.  Clusters  closer  than  this 
minimum  distance  are  combined  to  be  served  by  one  switch. 
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3.  LOCATING  RACK  HONK  SWITCHKP.  TN  A M1XK!)  Mill)  l A N.KTWORK 

In  tn is  section,  we  consider  a model  for  the  mixed  terrestrial 
satellite  network,  and  investigate  solution  approaches  to  the  corre 
spondincj  backbone  switch  location  problem. 

Ouj  network  model  is  as  follows: 

1.  The  terrestrial  network  consists  of  a set  of  store- 
and-forward  switches  (e.g.,  IMF's)  interconnected  by 
ground  channels  (a  distributed  subnet).  For  reliability, 
the  terrestrial  net  is  usually  required  to  be  2-connected. 

2.  The  ground  stations  are  colocated  with  the  satellite 
switches.  (This  assumption  is  used  here  ft  r convenience. 
However,  it  is  not  required  for  the  satellite  switch  se- 
lection procedure  developed  in  Section  <1 . ) Moreover, 

the  satellite  switches  form  a subset  of  the  set  of  store- 
and-forward  switches. 

3.  The  satellite  system  is  used  either  in  a dedicated 
access/broadcast  mode,  or  in  a multiple  access/broadcast 
mode.  Thus,  the  satellite  subnet  effectively  forms  a 
complete  graph, 

4.  All  the  switches  (regular  and  satellite)  have  the 
same  capacity  and  cost. 

5.  The  earth  station  cost  is  fixed,  while  the  satel- 
lite channel  cost  is  proportional  to  the  bandwidth. 

6.  The  low  delay  traffic  which  cannot  be  routed  over 
the  satellite  links  constitutes  only  a small  portion  of 
total  traffic  requirement,  and  has  higher  priority. 
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3.1  Backbone  Trans  in  i s s ion  C on  L Estimate  For  Mixed  Media  Networks 

Based  on  the  terrestrial  line  cost  estimate  developed  in  .Sec- 
tion 2.1 , we  can  approximate  the  backbone  network  cost  foj  a mixed 
media  network  as  follows:  For  any  two  switches  A,B  in  the  network, 

let 

= Cost  of  routing  a unit  flow  from  A to  B along 
terrestrial  links, 

- Cost  of  routing  a unit  flow  from  A to  B through 
some  satellite  links, 

= Nearest  satellite  switch  to  A.  ( S ( A ) may  coincide 
with  A. ) 

TFR 

Then  by  Eqs.  (1)  and  (2) , C^  ^ can  be  estimated  by 
TER 

C^  - C’T  x d(A,B)  + (estimated  unit  line  hardware  cost 
for  links  on  routes  from  A to  B) , (3) 

where  Cff  is  the  backbone  line  cost  approximation  factor  as  given 
in  Section  2.1  (in  fact,  C,p  = c.P/p,  where  c,P  and  p are  defined 
in  Section  2.1).  The  unit  line  hardware  cost  can  be  calculated 
by  estimating  the  number  of  links  on  an  average  path  between  A and 
B.  Since  the  hardware  cost  usually  constitutes  only  a small  frac- 
tion of  the  total  line  cost,  a very  rough  estimate  suffices. 

Specifically,  in  the  present  study,  the  hardware  cost  is 
estimated  in  the  following  fashion.  Similar  to  the  heuristic 
estimate  developed  in  Section  2.1,  for  a given  set  of  backbone 


TER 

LA,B 


,SAT 
A , B 


S(A) 
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switches  and  requirement  matrix,  we  estimate  the  channe 1 -m i 1 os 
and  the  number  of  channels  required  for  a terrestrial  backbone 
network.  From  this  wo  obtain  an  estimate  for  the  average  channel 
length,  L^ . The  number  of  links  on  an  average  path  between  two 
switches  A and  B is  then  estimated  by  V x d(A,B)/L  (I1  is  defined 
in  Section  2.1),  from  which  the  unit  line  hardware  cost  for  paths 
between  A and  B follows. 

The  unit  cost  of  routing  through  satellite  links  can  be  esti- 
mated by  (see  Figure  1) 


,.SAT  _ TER  TER 

"A,B  " LA,S(A)  B,S(B)  LS' 


(4) 


where  is  the  unit  satellite  bandwidth  cos'  . 

To  a first  degree  approximation,  for  the  low  priority  require- 
ment, the  route  with  the  lower  cost  is  preferred  (see  Figure  1),  and 
thus,  the  unit  cost  for  routing  from  A to  B,  ^ is  estimated  by 


CA,B  = min 


(C 


SAT 
A,  B' 


TER 
A , B 


(5) 


3 . 2 An  Approach  for  Locating  Backbone  Switches 


With  Eqs.  (3)-(5),  we  can  extend  the  MODULARIZED  ADD  algorithm 
to  select  the  satellite  and  terrestrial  switch  locations.  At  each 
iteration  of  the  ADD  algorithm,  we  not  only  can  select  one  more 
backbone  switch  location,  we  can  also  determine  whether  it  is  profit- 
able to  place  an  earth  station  at  that  location.  An  approach  to 
select  the  switch  is  as  follows:  First,  compare  the  (remaining) 

candidates  assuming  that  they  are  equipped  with  earth  stations.  Let 
the  best  candidate  chosen  be  and  corresponding  saving  S-^.  Then, 
compare  the  same  set  of  candidates  assuming  that  they  are  not  equipped 
with  earth  stations.  Let  the  best  candidate  chosen  be  C?2  and  corre- 
sponding saving  • The  CL  with  the  larger  saving  is  selected  as 
the  next  switch  location  (with  or  without  an  earth  station  at  the 
location  depending  on  whether  i is  1 or  2) . 
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S (A) 


S(B) 


SWITCHES  WITH 
EARTH  STATION 


SWITCHES  WITHOUT 
EARTH  STATION 


FIGURE  1:  ROUTING  ALTERNATIVES  IN  A MIXED  NETWORK 
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There  are  some  inherenL  weaknesses  with  this  approach,  aside 
from  computational  complexity  considerations.  For  example,  during 
some  iteration,  a switch  location  0 without  an  earth  station  may 
be  selected  on  the  basis  of  the  saving  criterion.  However,  when 
more  switches  are  selected  later,  it  may  become  more  profitable  to 
have  an  earth  station  at  Q.  Yet,  this  information  was  not  avail- 
able during  the  selection  of  Q. 

Alternatively,  one  can  use  some  very  simple,  but  fast  estimate 
for  the  backbone  line  cost.  The  following  is  an  outline  of  one 
such  approach.  Steps  1 and  2 are  used  to  estimate  the  incremental 
backbone  line  cost. 


1.  Rased  on  the  given  user  traffic  matrix  and  design 
constraints,  generate  a number  of  reasonably  "good" 
mixed  terrestrial -satellite  backbone  network  designs. 

2.  Based  on  the  designs  obtained  in  Step  1,  construct 
a function  which  approximates  the  backbone  lino  cost 
relative  to  the  number  of  backbone  nodes.  Let  this 
function  be  BC(N).  Also,  let 

ABC  (N ) = BC  (N+ 1 ) - BC  (N ) 

be  the  incremental  cost  function. 

3.  Use  the  ADD  algorithm  to  select  the  backbone  switch 

locations  with  the  following  modification:  Suppose  K 

switches  have  already  been  selected.  Then  the  saving 
that  can  be  achieved  by  selecting  a candidate  Q as  the 
(K-tl)-th  location  is 
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(local  iKtoss  saving  of  selecting  0 as 
t he  ( K+  1 ) - Lh  swi  t.ch  ; oca  t i on 

(fixed  cost  of  installing  a switch  at\ 

Q / 

- A13C  ( K ) . 

4.  Suppose  switch  locations  C, , C are  selected 

by  tne  modified  ADD  algorithm.  Determine  wnich  of  the 
C^'s  are  to  be  flipped  with  earth  stations  so  as  to 
optimize  the  backbone  network  cost. 

Based  on  experimental  results  presented  in  [NAC,  1976]  and  in 
Section  5.2,  we  found  that  for  terrestrial  networks: 

• The  number  and  location  of  backbone  switches  are 
rather  insensitive  to  variations  in  unit  backbone 
cost . 

• The  best  selection  made  by  the  basic  ADD  algorithm 
given  in  Section  2.2  is  only  slightly  better  than 
the  best  selection  made  by  the  simplified  scheme 
described  in  Steps  1-3. 

We  thus  conclude  that  for  a general  backbone  network  (terrestrial, 
mixed  terrestrial-satellite,  hybrid  packet-circuit,  etc.),  the  pro- 
cedure based  on  Steps  1-3  is  sufficient  to  obtain  a good  selection 
of  backbone  switches. 

Step  4,  the  optimization  of  number  and  location  of  satellite 
switches,  is  in  itself  a rather  difficult  problem.  It  is  discussed 
in  the  next  section. 


8k(Q) 
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4.  SATELLITE  SWITCH  SELECTION 


In  this  section,  wo  consider  the  problem  of  selecting  the 
satellite  switch  (and  thus  tin-  associated  ground  station)  locations 
from  among  the  backbone  switch  locations.  More  formally,  the  prob- 
lem is  as  follows: 


(h_vqn: 

• Switch  locations. 

• Switch- to-switch  traffic  requirement. 

• Ground  station  cost  (including  the  cost  of  con- 
nection to  the  terrestrial  network). 

• Satellite  bandwidth  cost. 

• Satellite  access  technique. 

• Terrestrial  bandwidth  cost. 

Optimize : 

Total  communication  cost  D: 

D = ground  station  cost  + terrestrial  trunk 
cost  + satellite  bandwidth  cost. 
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Over  the  Vnr  i ib]  es  : 

• Number  and  location  of  satellite  switch.es  (ground 
stations ) . 

• Terrestrial  network  topology. 


Such  that: 


• Traffic  requirements  are  met. 

• Appropriate  constraints  (delay,  reliability,  etc.) 
are  satisfied. 

• Switch  capacity  constraints  arc  met. 

To  simplify  the  solution  procedure,  we  use  the  following  as- 
sumptions; they  can  be  relaxed  by  modifying  the  basic  procedure 
appropriately : 

• Switch  capacity  is  unlimited.  This  is  quite  realistic 
if  we  use  modularized  switches  such  as  the  Pluribus 
IMP  (HEART,  1973]. 

• Switch  cost  varies  linearly  with  capacity,  and  the 
fixed  cost  portion  dominates.  Again,  this  is  a 
reasonable  cost  model  for  the  modular  type  switch. 

• The  traffic  J assumed  to  consist  of  two  priorities; 
the  average  delay  for  the  high  priority  traffic  is 
smaller  than  the  propagation  time  over  the  satellite 
channel  and  hence  must  be  accommodated  by  the  ter- 
restrial subnet.  The  low  priority  traffic  can  tolerate 
satellite  propagation  delay. 
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In  the  solution  procedure,  the  delay  of  the  low  priority 
tralfic  is  not  cxplicity  taken  into  account ; it  in  implicitly 
used  in  the  average  link  utilisation  in  deriving  the  cost  of 
routing  a unit  flow  over  a channel. 

4.1  B asi  c:  Soli  1 1 ion  A| >proach 

From  the  transmission  cost  cstima4 os  developed  in  the  last 
section  (Kqs . ( 3 ) - ( 5 ) ) , we  can  approximate  the  backbone  cost  for  any 
switch  set  configuration.  This  forms  the  basis  of  our  exhaustive 
search  algorithm. 

TEH 

First,  notice  that  ^ , the  line  cost  of  sending  a unit  flow 

from  A to  B,  can  be  regarded  as  a constant  independent  of  the  satel- 
lite switches  selected. 

Now  let  l be  a subset  of  the  switches.  Suppose  the  nodes  in  f 

are  used  as  satellite  switches.  Then,  the  satellite  transmission 

SAT 

cost  for  a unit  flow  between  A and  B,  C^‘^(4>),  can  be  calculated 
from  Eq . ( 4 ) . The  unit  transmission  cost,  C ($»),  can  then  be 

f o 

calculated  from  Eq . (5).  It  fellows  that  the  total  backbone  com- 
munication cost  (excluding  the  switch  cost)  for  the  low  priority 
traffic  is: 


**  I Xv'GRD  + 


where 


= Ground  station  cost, 

vjKL/ 

t^  ^ = Low  priority  requirement  from  A to  B. 

The  optimum  satellite  switch  set  can  thus  be  determined  by 
selecting  the  collection  !>  with  minimum  cost  Cbf). 
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For  a given  sot  of  N switcher,,  the  total  number  of  poss.ih.Io 
N 

suoscts  is  2 . Also,  for  each  subset  T,  the  computation  of  C(!) 

requires  on  the  order  of  operations.  Thus,  the  complexity  of  title 

?.  N 

above  procedure  if;  on  the*  order  of  N‘  x 2 . Since  for  most  practi- 
cal problems,  N is  small  (on  the  order  of  10),  thus  the*  computation 
of  the  optimum  subset  li’  is  feasible.  Moreover,  as  will  be  seen  in 
Sections  4.2  and  4.3,  there  are  various  ways  to  reduce  the  computa- 
tion. 


4 . 2 Cyclic  Sub  si » t_  Order  i iig 


Suppose 
id  f 3 

(4*)  }, 


and 


are 


tained  from  4>  by  deleting 

^ C A , B ^ 1 A , 13 ' and  f(A,U<:) 
for  T are  in  storage. 


t*wo  sets  of  switches  such  that  (j>'  is  ob- 
a switch.  Wo  investigate  how  to  compute 
A , given  that  the  corresponding  items 

/\  f ll 


C ( t) 


First,  we  assume  that  for  any  two  switches  A and  13  in  a net- 
work, the  number  of  links  (and  lienee  the  line  hardware  cost)  between 
A and  B is  proportional  to  the  distance  between  A and  B.  It  follows 
that  for  any  four  switches  A,  B,  C,  and  D 


d (A, B)  < d (C, D) 


implies 


,TER  TER 
A , 13  1 LC,D 


(7) 


TER 

Moreover,  the  triangular  inequality  holds  for  C , i.e.,  for  any 
three  switches  A,  B,  and  C, 


TER  , fER  TER 

A, B LB,C  - LA,C 


(0) 


Let  II  (Q)  be  the  set  of  switches  that  has  Q as  the  nearest 
satellite  switch  in  t . Given  two  switches  A and  B,  note  that 
CA  ')  could  be  different  from  C’A  only  if  either  A or  B is 

in  Ii  (Q)  . Suppose  A is  in  ii  (Q)  . Let  S'  (A)  be  the  satellite  switch 
in  <J>'  nearest  to  A.  Note  that 


d (A, S ' (A) ) > d(A,Q) 


H6nce,  it  follows  from  Eqs.  (7)  and  (4), 
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rS AT  , 
CA,I3(i  } 


• cSAT  ( <:> 

A,Bl 


We  conclude  that; 


1. 


lf  CA,B(1>)  “ c™.  then 


CA , 11 ( 41 ' * CA  , B CA,B(t' 


(Note  that  if  A , B c.  II  (Q)  , then  C (<I>) 

l\  § 13 


^TERv 

l'A,B  ' 


2.  lf  S' (A)  = S' (B) , then 


TKR 

CA,B<$,)  = CA,B 


where  S'(B)  is  the  nearest  satellite  switch  to  B in  , 


3.  Otherwise, 


cSATK')  = cTKR  + cTE“  + c 

A,BV  ' A, S' (A)  B,S’(B)  S 


and 


„ , . . r „TER  „SAT  , , , x , 

CA / B ^ * mln  CA , B ' CA , B ^ 

From  the  discussion  in  the  last  paragraph,  C ( 1 ' ) can  be  com- 
puted from  C(fl>)  as  follows: 

1.  Set  C(4>')  = C($)  . 

2.  Fo^  each  A c II  (Q)  , calculate  S'  (A). 


SAT 

3.  For  each  A c II  (Q)  and  each  B,  compute 
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4.  For  each  A * II  (Q)  and  each  switch  B / 11  (Q)  , if 


CA,D('1')  ^ Ca'/u'  thon 


a . Compute  ^ (<I  1 ) 


b.  C(<1'')  = C (<J>  ’ ) CA>IJ  (<!'')  - CA>Ij(<M 


5.  C(*')  = C<*')  - CGRD. 

Next,  we  consider  the  situation  when  1>  ' is  obtained  from  f by 

adding  a switch  Q not  in  4» . Let  n (Q ) be  the  set  of  switches  that 

has  Q as  the  nearest  satellite  switch  in  v ' . Similar  to  before, 

C (4* ' ) v;ould  be  different  from  C (T)  only  if  either  A or  B is 
i\  f 13  /\  z 1J 

in  11  (Q)  . For  each  switch  If,  let  S (B)  be  the  satellite  switch  near- 
est to  B in  4> . Let  A be  a switch  in  ..  (Q)  . Then  we  have 


d (A,Q)  ' d (A , S (A ) ) , 


and  hence  by  Eqs . (7)  and  (4), 


CSAT(*  ' ) < CSAT  (']>)/ 
A,B^  ' - ^A,BU  ' ' 


for  any  switch  B.  Note  that  if  B is  also  in  II  (Q)  , then  by  Eq.  (8)  , 


TEK  TER  TER  TER  TER 

^A,B  - UA,Q  U ,Q  - ^A , S (A ) B,S(B) 


Hence , 


CAfB(*' 


= CA,B(*> 


rTER 

A,B’ 
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We  thus  conclude  that 


i K 

CA,B(i' 

) /■ 

CA,B(1,)  °,,ly  K 

B / il  (Q)  . 

1 2- 

lf  CA,B 

(*) 

= then 

1 H 

CA.B(** 

) ' 

SAT  / , , v , 

ca  , n ' ’ * ' and 

1 ■ 

CA,B<*' 

) “* 

C (S’)  - C1LR 

^A , B A,Q 

TER 

A,S(A) ’ 

3. 

Otherwise , 

1 III 

CSAT. 
*"A , B ' 

) = 

TEH  TER 

A , Q C,S(B) 

+ cs. 

^ and 

LA,Biq 

) = 

mm  {. C A # I5  ( ^ ) * 

..SAT  , , , , 

CA,B("’ 

From  the  discussion  in  the  last  paragraph,  C(<I>')  can  be  com- 
puted from  C($)  as  follows: 

* 

1.  Set  C ($  ' ) = C(4»)  . 


2.  Find  the  set  n (Q)  of  switches  that  have  Q as  the 
nearest  satellite  switch  in  <\> ' . 

SAT 

3.  For  each  A c II  (Q)  and  each  B,  compute  (<P ' ) . 

4.  For  each  A c II  (Q)  and  B / II  (Q)  , 

SAT 

a.  If  CA#B(*)  - CA#B(*},  then 
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1 . 

C (cj  ’ 

A,B  v 

SAT 
A , B 

(4  1 ) , 

2. 

C (4-  1 ) 

G(t 

’)  + 

TER 

A,Q 

TER 
A,S  (A) 

Otherwi sc , 

1. 

C (4>  ’ 

A,Biq 

) - 

min 

O' 

c Arr 

■),  c™>-) 

2. 

c(«r>  = 

C (T 

’)  + 

cA,n>* 

’>  - CA,U<;> 

5.  C ( T * ) ~ C ( 4»  ' ) 4 CGRD. 

From  above  discussions,  it  follows  that  substantial  savings  in 
cost  compulation  can  be  obtained  if  each  successive  set  of  satellite 
switches  is  obtained  irom  the  previous  set  by  either  adding  a switch 
or  deleting  a switch.  Suppose  there  are  N switches  in  the  backbone 
network.  Then  each  subset  S of  switches  corresponds  to  a unique  li- 
bit binary  word  9nc*n-1'  •••*  where  g^  is  1 if  and  only  if  switch 

i is  in  S.  Hence,  what  we  desire  is  an  ordering  of  the  N-bit  words 
such  that  any  two  successive  words  differ  in  one  and  only  one  bit. 

A binary  code  with  such  a property  is  called  the  Gray  code  or  the 
cyclic  code  [BOOTHROYD,  1964  ],  f GSCHWTND , 1975].  A discussion  on 
the  generation  procedure  of  one  type  of  such  code  words,  the  re- 
flected binary  cyclic  code,  is  presented  in  Appendix  A. 

4 . 3 Complexity  of  the  Cycl ic  Subset  Ordering  A 1 gor d Lhm 

Below,  we  estimate  the  time  complexity  of  the  procedure  pre- 
sented in  Section  4.2.  Suppose  there  are  k satellite  switches. 
Then,  on  the  average,  each  satellite  switch  is  the  nearest  satellite 
switch  for  N/k  switches,  where  N is  the  total  number  of  switches. 
Hence,  the  total  time  required  over  all  the  subsets  of  N switches, 

'"Total ' 1S 
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tf  (Time  required  to  generate  ail  subsets) 

1 OtO  1 

N 

+ Y,  [(number  of  k-element  subsets) 
k 0 


x (average  time  required  lor  each  k subset)] 


= (time  required  to  generate  all  subsets) 


+ 


+ 


N 

L 

k=l 


U)  . 


From  Appendix  A,  the  generation  of  all  subsets  such  that  each 

two  successive  subsets  differ  in  one  and  only  one  element  can  be 

"M*f  1 

done  in  time  0(2  ).  Hence/ 


tTotal 


- 0(2N+1) 


+ N 


H 

>: 

k-  3. 


<i* 


N 

k 


N. 


= 0(2 


N-!  1 , 


2 N 
+ N + N £ 

k=l 


<;> 


x 


N 

k * 


We  will  show  that 


N 

l 

k=l 


■:> 


3 


x 


t 


and  consequently/ 

N-f  1 2 N 

Wax  = 0(2  > + N + 3N  x 2N 

= 0 (N  x 2N) 

This  is  an  o^der  N improvement  over  the  straight-forward  ap- 
proach presented  in  Section  4.1. 
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ASSERTION 


v N 0N 

>-  (k)  x k *'  3 x 2 ' 

1 K k “ 


PROOF 


N M ,, 

v / N x W 

1 (V)  x k 

k=l  K k 


N 

>:  ( 

k=l 


N+l  k+1  N 
k+1  X N + l X k 


N 


N+l, 


k=l 


k+1 


N 


k + 1 


N+l 


v /N+l.  N #1,1. 
k+1  X N+f  A (1  k} 

k~l 

= x nTTx  [(14>  + <1+N=k)] 


T’o,  « 

k=l  K 1 


N+l  * (2+|fN-k>- 


Nov;,  notice  that  for  1 £ k < [N/2]  , the  maximum  value  for  ~ + — ^ 
occurs  at  k-1.  Thus 


_N_ 

N+l 


x 


(2+±+-l_)  < JL. 

U k N-k'  - N+l 


X (2  + l+^__^) 


N 3N-2 
N+l  X N-l 


< 3. 


But,  then  we  obtain 
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N 

(N/2) 

N 

X k 

>: 

k-  1 

(ku’ 

X INH1 

[N/2] 

V 

k=  1 

(N+1) 
k -f  .1  ' 

x 3 

< 

3 x \ 

0N+ 
x 2 

1 

N 

=3x2.  Q.E.D. 

4 . 4 Other  Problem  Reduction  Techniques 

The  procedure  described  in  Section  4.2,  although  is  an  improve- 
ment over  the  straightforward  procedure,  still  takes  time  N x 2N. 
Hence,  for  switch  sets  of  moderate  to  large  size  (say  20  to  50 
switches),  other  problem  reduction  techniques  are  needed  to  obtain 
a solid  ion.  In  the  following,  we  describe  several  such  possibili- 
ties . 

First,  with  the  given  earth  station  cost,  terrestrial  trans- 
mission cost,  traffic  matrix,  etc.,  we  can  obtain  various  upper  and 
lower  bounds  on  the  number  of  satellite  switches  that  can  bo  de- 
ployed in  a region.  For  example,  the  following  result  is  derived 
in  Appendix  B. 


ASSERTION 

Suppose  A is  selected  as  a satellite  switch.  Then,  for  any 
switch  B with 

t x CTER  < C 
B A,B  - uRD' 

it  is  not  cost-effective  to  implement  a satellite  switch  at  B, 
(where  t^  is  the  total  low  priority  traffic  requirement  at  B) . 

These  types  of  bounds  are  especially  useful  when  the  number  of 
switch  sites  is  large. 
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Second,  during  experiments,  we  found  that  a number  of  the  loca- 
tions consistently  qet  picked  as  satellite  switch  sites,  independent 
of . the  total  number  of  backbone  switches,  and  within  a wide  range  of 
cost  and  traffic  variations  (o.q.,  TTNKPK,  McCKLLAN,  ALBANY).  These 
are  usually  nodes  with  high  traffic  requirements.  Hence,  during  the 
satellite  switch  selection  process,  we  can  include  these  nodes  in 
the  satellite  switch  list,  and  exclude  some  of  their  neighboring 
nodes  from  consideration. 

Moreover,  from  the  experimental  results  reported  in  Section  5.6.1, 
we  observe  that  for  a given  Host  and  terminal  data  base,  the  number 
of  satellite  switches  selected  is  almost  independent  of  the  size  of 
the  backbone  switch  set.  It  follows  that  for  a design  involving  a 
large  number  of  switches,  we  can  first  estimate  the  number  of  satel- 
lite switches  by  applying  the  procedure  to  a problem  with  much  fewer 
switches,  and  then  limit  the  search  to  only  switch  sets  with  size 
close  to  the  estimated  value. 
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5 . EXPER  1 MENTAL  RESULTS 

5.1  Sy s t em_  Mod e I and  Raso  P a r a me t o r Values 

Wo  have  conducted  extensive  experiments  with  both  t.he  satellite 
switch  selection  algorithm  developed  in  Section  4,  and  the  MODULAR- 
IZED ADD  algorithm,  on  the  AUTODIN  II  system  with  the  AUTOVON  sites 
as  candidate  backbone  switch  locations.  The  AUTODIN  11  and  AUTOVON 
system  size  is  as  follows: 

AUTODIN  II  SYSTEM 


Total  Number  of 

User  Locations 

- 

300 

Number  of 

Host  Computers 

= 

86 

Number  of 

Conc'-ntrators 

■- 

26 

Number  of 

TDEX ’ s 

= 

00 

ei 

Number  of 

TCU  ’ s 

- 

11 

Number  of 

Isolated  Terminals 

= 

101 

Total  Traffic 

= 

1 

AUTOVON  SYSTEM 

Total  Number  of 

AUTOVON  Sites 



60 

Megabits/Sec . 


Figures  2 and  3 are  the  AUTODIN  II  and  AUTOVON  locations,  respec- 
tively. 

The  assumptions  used  in  the  mixed  media  backbone  network  design 
are  as  follows: 


Proportion  of  Low  Delay  Traffic 
Average  End-to-End  Delay  For  Low 
Delay  Traffic 
Protocol  Overhead 


= 10?, 

100  msec. 
= 4 5% 
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FIGURE  2:  USER  SITES  IN  THE  AUTODIN  II  SYSTEMS 
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Routing  Over  v1*  id  = 11 

Average  I’.iduT  Length  - 500  bits 

Routing  Method  = Optimal  MuJ tipath 

Satellite  Channel  Access  Technique  - Dedicated  Up-Link/ 

Broadcast  Down-Link 

Tables  1 and  2 list  cost  and  characteristics  of  local  access 
and  backbone  facilities  used  in  the  design. 

5 . 2 Sen  si  t ivi  ty  Stu  d ios  on  the  Backbone  I ,o  c a_t  ion  of  Terr  o s t r i a 1 
Not  works 

Here,  we  present  some  experimental  results  on  the  backbone 
switch  location  problem  for  terrestrial  networks . The  main  purpose 
of  this  section  is  to  provide',  in  some  sense,  experimental  justifi- 
cations for  the  switch  location  approach  given  in  Section  3.2. 

We  compare  the  following  two  different  modifications  of  the 
ADD  algorithms  on  terrestrial  network  backbone  designs: 

1.  The  procedure  presented  in  Section  2.2  (and  [NAC, 

1976])  which  uses  detailed  backbone  line  cost  tradeoff 
evaluations.  We  call  this  the  DB  (Detailed  Backbone 
Estimation)  approach. 

► 

2.  The  procedure  presented  in  Section  3.2  which  uses 
only  rough  incremental  backbone  line  cost  estimate, 

ABC(K),  dependent  on  the  number  of  backbone  switches 
used.  We  call  this  the  RB  (Rough  Backbone  Estimation) 
approach. 

For  simplicity,  in  this  study  wo  make  ABC(K)  the  same  for  all 
values  of  K.  Thus,  we  use  only  ABC  to  denote  ABC(K).  From  exper- 
ience with  the  AUTODIN  II  system  design,  the  backbone  line  cost  varies 
approximately  linearly  with  the  number  of-  backbone  switches  within 
the  range  of  5-20  switches,  [NAC,  1976],  thus  the  above  approxima- 
tion is  justified. 
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CIL1TY  TYPE 

COSTS  (3 /MO) . 

KB/S  L i lie 

j 

2 . 0/mi ] e 

KB/S  Line 

2 . 0/rni lc 

KB/S  Line 

2 . 0/mi Ic 

KB/S  Line 

12 . 0/mi lc 

KB/S  Line 

100. 

KB/S  Line 

240. 

KB/S  Line 

570. 

KB/S  Line 

CD 

tn 

o 

Isolated  Terminal 

TDMX 

TCU 

Concentrator 


TDMX 

TCU 

Concentrator 


0 

67, 

1170 
272  4. 


13, 

0 

0 


TABLE  1:  LOCAL  ACCESS  FACILITY  UNIT  COST 
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FACILITY  TYPE 

COST  ( $/MO . ) 

Software  Development  One-Time  Fixed  Cost 

60,000 

Satellite  Earth  Station  Cost 

10,000 

Satellite  Channel  (50  KB/S-Simplex) 
Space  Segment  Cost 

500 

Satellite  Channel  (50  KB/S-Simplex) 
Access  Cost 

300 

Terrestrial  Channel  (50  KB/S] 
Mileaqe  Cost 


5/ mi le 


Terrestrial 

Termination 

Channel  (50 
Cost 

KB/S) 

Switch  Node 

(Pluribus ) 

Fixed  Cost 

Trunk  Line 

Host  or  Concentrator  Line 
TCU  Line 
MUX  Line 

2.4  KB/S  Terminal  Line 
4.8  KB/S  Terminal  Line 
9.6  KB/S  Terminal  Line 
50.  KB/S  Terminal  Line 


425/end 
10 , 600 


TABLE  2:  BACKBONE  FACILITY  UNIT  COST 
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Table  3 lists  the  cost  estimates  of  the  designs  generated  by 
the  two  procedures.  (Although  each  set  of  switch  locations  is 
generated  by  using  a particular  ABC  value,  in  order  to  have  an  adequate 
comparison,  the  backbone  costs  listed  in  Table  3 are  estimated  by  the 
same  formula,  Eq.  (2),  with  unit  costs  given  by  Table  2. 

Comparing  the  design  costs  in  Table  3;  we  can  conclude  that 

1.  The  cost  of  the  best  (lowest  cost)  designs  for  dif- 
ferent values  of  ABC's  differ  very  little  from  one  an- 
other. This  implies  that  the  exact  determination  of  ABC 
is  not  important.  A rough  estimate  for  ABC  can  usually 
give  us  a reasonably  good  selection  of  backbone  switches. 

2.  The  best  designs  obtained  with  the  RB  approach  are 
slightly  worse  than  the  best  design  obtained  with  the  DB 
approach  (the  difference  is  less  than  1.1%).  This  sug- 
gests that  a coarser  approach  such  as  RB  suffices  for 
most  applications.  The  RB  approach  is  inherently  much 
faster  than  the  DB  approach. 


In  Table  4,  we  compare  the  designs  listed  in  Table  3 that 
have  the  same  number  of  backbone  switches  but  with  different  total 
cost.  We  can  observe  the  following: 


1.  For  the  designs  giving  the  same  number  of  backbone 
nodes,  the  design  with  lower  local  access  cost  almost 
always  also  has  lower  total  cost. 


2.  For  the  designs  giving  the  same  number  of  backbone 
nodes,  the  design  with  lower  local  access  cost  usually 
has  higher  backbone  line  cost. 
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DESIGN  APPROACH 

RB 

DB 

Number  ABC 

of  Backbone  (K$/mo) 

Nodes 

0.0 

6.6 

12.0 

-- 

7 

757 

757 

757 

75/ 

8 

- 

- 

*745 

753 

9 

*749 

751 

111 

750 

10 

- 

- 

- 

11 

- 

*750 

750 

*742 

12 

757 

757 

761 

756 

13 

- 

- 

- 

- 

14 

784 

- 

- 

791 

13 

- 

- 

830 

- 

16 

- 

844 

- 

797 

17 

- 

- 

“ 

- 

18 

858 

— 

~ 

828 

TABLE  3:  ESTIMATED  TOTAL  COST  OF  DESIGNS  GEN  I)  RAT  ED 

BACKBONE  NETWORK  COST  ESTIMATED  BY  EQUATION  (2) 
UNIT  = K$/mo 


^Lowest  cost  designs 
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NUMBER  OF 
BACKBONE 
• NODES 

SOLUTION 

PROCEDURE 

A BC 

( K $ / MO ) 

LOCAL  ACCESS 
COST  ( K$/MO , 

BACKBONE  LINE 
COST  (K$/MO) 

TOTAL  COST  j 
( K $ / MO ) 

8 

DB 

- 

429 

164 

753 

RI3 

12.0 

416 

169 

745 

9 

DB 

- 

406 

173 

750 

RB 

0.0 

383 

184 

749 

RB 

6.6 

406 

173 

751 

RB 

12.0 

417 

182 

771 

11 

DB 

- 

360 

190 

742 

RB 

6.6 

365 

192 

750 

RB 

12.0 

365 

192 

750 

12 

DB 

- 

357 

1 96 

756 

RB 

0.0 

356 

198 

757 

RB 

6.6 

356 

198 

757 

RB 

12.0 

360 

197 

761 

14 

DB 

- 

362 

205 

791 

RB 

0.  0 

348 

211 

784 

■ 

15 

DB 

- 

329 

223 

797 

RB 

379 

218 

830 

18 

DB 

- 

326 

235 

828 

RB 

0.0 

354 

237 

858 

TABLE  4:  COMPARISON  OF  DESIGNS  WITH  SAME  NUMBER  OF  SWITCHES 

(BASED  ON  DESIGNS  LISTED  IN  TABLE  3) 

BACKBONE  NETWORK  COST  ESTIMATED  BY  EQUATION  (2) 
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From  the  above  observations,  we  conclude  that  for  a given  num- 
ber of  backbone  switches,  a cjood  hcuiistic  cri  I rnon  in  selecting 
low  cost  designs  i : to  minimize  the  local  access  cost.  This  agrees 

with  the  intuition  that,  because  the  .local  access  cost  as  the  larg- 
est component  of  total  cost,  it  has  the  greatest  influence  on  the 
total  design  cost. 


5 . 3 Accuracy  of  Backbone  Network  Cos t Estimates 

The  ADD  algorithm,  given  in  Section  3.2,  for  selecting  the  back- 
bone switches  (for  mixed  networks)  estimates  the  backbone  incremental 
cost  of  one  additional  switch,  given  that  there;  are  already  K switches 
selected,  by  a single  value:  ABC(K).  The  algorithm  tjiven  in  Section  4 

for  selecting  the  satellite  switches  estimates  the  backbone  network 
cost  by  Lq . ( 6 ) : 


cuo  = |*|  X cGRD  + tA(B  X cA|B(*), 


where 


$ = the  set  of  satellite  switches, 


C«n~  = ground  station  cost, 
wKU 


t = low  priority  requirement  from  A to  13, 
f 15 


and 


CA  B ( 4> ) is  given  by  Eqs.  (3)  - (5). 

Here  we  examine  the  accuracy  and  the  consistency  of  these  backbone 
network  cost  estimates. 
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We  have  manually  designed  several  m j xed~med  i a networks  with 
varying  numbers  of  backbone  nod  or; . Based  on  these  designs,  wo 
estimate  the  backbone  network  incremental  cost  (excludincj  the 
switch  cost)  to  be 

ABC ( N ) $4 ,500/month,  (9) 

and  the  backbone  network  cost  (excludincj  the  switch  cost)  to  be 

BC (N)  $88, 100/month  + ABC(N)  x N,  (10) 

for  5 <_  N £ 20.  Several  sets  of  backbone  switch  locations  were 
then  selected  using  the  ADD  algorithm  based  on  these  values  over 
a range  of  alternatives.  For  each  set  of  the  backbone  switches 
selected,  the  optimum  satellite  switch  locations  were  then  selected 
using  the  satellite  switch  selection  algorithm  based  on  Eg.  (G). 

Actual  topological  network  designs  were  then  generated  using  these 
satellite  and  terrestrial  switch  locations.  In  Figure  4,  we 
compare  the  backbone  network  cost  of  the  actual  mixed  network 
designs  to  the  estimates  obtained  from  Eg.  (10)  .‘  As  can  be  seen, 
the  estimates  and  the  actual  costs  are  uniformly  close  (always  with- 
in 10%) . The  fact  that  the  optimum  number  of  satellite  switches  is 
fairly  independent  of  the  switch  set  size  (see  Section  5.G.1)  may 
contribute  partially  to  the  close  fit  of  the  estimates.  It  is  some- 
what surprising  that  the  backbone  network  cost  (as  a function  of  the 
number  of  switches)  can  be  approximated  so  closely  by  a straight  line. 
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BACKBONE  NETWORK 
COST  (K?/MO) 


200 


BACKBONE  COST  ESTIMATE  USED  IN  ADD  (BC(N)) 


BACKBONE  COST  CURVE  ESTIMATED 
BY  ACTUAL  NEAR  OPTIMUM  DESIGNS 


50 


O ESTIMATED  OPTIMUM  MIXED  NETWORK 
BACKBONE  COST 

A -ACTUAL  DESIGN  NETWORK 
BACKBONE  COST 


VvV-, | j j p 

6 7 8 9 10 


11 


12 


I 

13 


I 

14 


I 

15 


I 

16 


NUMBER  OF 
SWITCHES 


I 

17 


► 


FIGURE  4:  COMPARISON  OF  ACTUAL  BACKBONE  NETWORK  COSTS 

(EXCLUDING  SWITCH  COST)  VS,  ESTIMATES  DY 
EQUATION  (10)  USED  IN  ADD  ALGORITHM 
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In  order  to  assess  the  accuracy  of  Eq . (G)  and  the  satellite 

switch  selection  algorithm,  we  consider  an  11 -switch  set  selected 
by'  the  ADD  algorithm.  For  a given  k,  O^k  11,  wc  use  the  algorithm 
to  select  the  k optimum  satellite  switches.  A mixed  network 
topological  design  is  then  generated  based  on  the  11  switches  and 
these  selected  satellite  switches.  In  Figure  5 we  compare  the 
backbone  network  cost  of  the  actual  design  with  a given  > : of 
satellite  switches  to  the  estimates  by  Eq.  (6)  on  the  same  set  of 
satellite  switches.  It  can  be  seen  that  there  is  a gap  between 
the  actual  costs  and  the  estimates  (about  10-15%  difference).  How- 
ever, the  difference  is  quite  consistent,  and  the  general  pattern 
of  the  two  cost  curves  arc  similar.  We  thus  conclude  that  10 q . (G) 

is  adequate  for  the  purpose  of  satellite  switch  site  selection. 

5.4  Total  Network  Cost  Versus  Number  of  Backbone  Nodes 


Figure  6 plots  the  total  network  cost  estimates  (in  which  back- 
bone costs  are  estimated  by  Eq.  (10))  as  a function  of  the  number  of 
backbone  switches.  Observe  that,  similar  to  terrestrial  network 
designs,  these  points  represent  very  closely  a convex  (U-shaped)  curve, 
and  this  curve  is  quite  flat  near  the  optimum  point.  Moreover,  the 
mixed  networks  usually  offer  considerable  savings  over  the  terrestrial 
networks . 

For  illustrative  purpose,  the  following  designs  are  plotted 
in  Figures  7a-7f: 

a.  A terrestrial  network  design  for  8 prespecified 
backbone  switches,  with  a total  system  cost  of  $764K/mo. 

b.  A mixed  network  design  for  the  same  8 prespecified 
backbone  switches,  with  a total  system  cost  of  $722K/mo. 

c.  A terrestrial  network  design  for  11  ADD  algorithm 
selected  switch  locations  with  a total  system  cost  of 
$734K/mo. 
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BACK MONK  COST 
(K$/MO) 

A 


100 

□ ESTIMATED  OPTIMAL  DESIGN  COST  BY 
EQUATION  (G) 

A ACTUAL  DESIGN  COST 

50  - 


1 


2 


T 1 T 


NUMBER  OF 
SATELLITE 

SWITCHES 

* 


3 4 5 6 


8 9 10  11 


FIGURE  5:  COMPARISON  OF  BACKBONE  NETWORK  COSTS  (EXCLUDING  SWITCH 

COST)  VS.  ESTIMATES  BY  EQUATION  (0)  USED  IN  SATELLITE 
SWITCH  SELECTION  ALGORITHM 
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TERRESTRIAL  BACK SON 


UMERALS  ALONG  THE  LINKS  INDICATE  THE  NUMBER  OF  50KB/S  CHANNELS  USED 
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THE  NUMERALS  ALONG  THE  LINKS  INDICATE  THE  NUMBER  OF  50K3/S  CHANNELS  USED 
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NOTE:  THE  NUMERALS  ALONG  THE  LINKS  INDICATE  THE  NUMBER  OF  50KB/s  CHANNELS  USED 
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d.  A mixed  network  design  for  the  same  11  ADD  algo- 
rithm selected  switch  locations  with  a total  system  cost 
of  $607K/mo. 

e.  A terrestrial  network  design  for  16  ADD  algorithm 
selected  switch  locations  with  a total  system  cost  of 
$764K/mo . 

f.  A mixed  network  design  for  the  same  16  7\DD  algorithm 
selected  switch  locations  with  a total  system  cost  of 
$726K/mo . 


Observe  that  the  mixed  network  topoloaies  is  usually  much  simpler 
than  the  corresponding  terrestrial  network  topologies.  The  terres- 
trial portion  of  a mixed  network  usually  consists  of  a low  capacity 
loop  through  all  the  switches  plus  a few  more  cross  links. 

5.5  Backbone  Network  Cost  Versus  Number  of  Satellite  Switches 


Figure  8 plots  the  estimated  backbone  network  cost  (cost 
estimate  based  on  Eq.(7))  of  the  best  designs  with  various  given 
number  of  satellite  switches,  with  different  backbone  switch  sets 
and  throughput  levels.  As  can  be  seen,  the  total  cost  curve  is 
generally  very  flat.  This  fact  is  of  great  help  in  making  near- 
optimal  designs. 

In  Figure  9,  the  five  backbone  component  costs  (earth  sta- 
tion cost,  satellite  space  segment  cost,  satellite  access  cost, 
terrestrial  mileage  cost  and  terrestrial  hardware  cost)  are  plot- 
ted as  functions  of  the  number  of  satellite  switches  for  the  designs 
generated  with  the  11  backbone  switch  set.  It  can  be  seen  that  as 
more  satellite  switches  are  used,  the  cost  of  the  terrestrial  com- 
ponents decreases,  while  the  cost  of  the  satellite  components  in- 
creases. This  is  because  aside  from  the  earth  station  cost,  satel- 
lite transmission  is  usually  much  cheaper  than  terrestrial  trans- 
mission. Hence,  as  more  satellite  switches  become  available,  more 
and  more  terrestrial  channels  are  replaced  by  satellite  channels. 
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TERRESTRIAL  HARDWARE 
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K EARTH  STATION 
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TERRESTRIAL  MILEAGE 


SATELLITE  SPACE 
SEGMENT  COST , 


fV-. 

-0 0- 


^ SATELLITE  ACCESS  COST 


NUMBER  OF 
SATELLITE 
SWITCHES 


FIGURE  9 : BACKBONE  COMPONENT  COST  VS.  NU MBER  OF  S ATE ELITE 

SWITCHES , 11  BACKBONE  SWITCH  SET , 100%  THROUGHPUT 


NOTE:  Network  Cost  Estimated  by  Eq.  (6) 
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Moreover,  as  is  obvious  from  Figure  9,  the  earth  station  cost 
and  the  terrestrial  mileage  cost  are  the  two  main  cost  components, 
and  compared  with  the  other  cost  components,  they  are  much  more 
sensitive  to  the  variation  in  the  number  of  satellite  switciies. 
Thus,  the  main  tradeoff  in  determining  t he  number  of  satelli t e 
switches  is  between  the  earth  station  co st  and  the  Leriusf r i al 
mileage  cost . 

5.6  Sensitivity  Study  on  the  Optimum  Number  of  Satellite  Switches 
5.6.1  Effect  of  Backbone  Switch  Set  Size 


On  the  switch  sets  we  have  investigated  (all  except  the  8-node 
switch  sets  are  obtained  by  the  ADD  algorithm),  the  optimal  number 
of  satellite  switches  are  as  follows: 

Backbone  Switch  Set  Size  Satellite  Switch  Set  Size 


7 4 

8 5 

9 4 

11  4 

13  4 

16  4 


It  thus  appears  that  the  number  of  satellite  switches  required  for 
a good  design  is  independent  of  the  total  number  of  switches,  given 
that  the  switch  sites  are  chosen  optimally,  e.g.,  by  the  ADD  algorithm. 
This,  together  with  the  fact  that  the  total  backbone  network  cost  curve 
is  very  flat  ' vith  respect  to  the  number  of  satellite  switches), 
suggests  that  in  selecting  the  satellite  switches,  one  needs  only 
consider  switch  sets  with  size  within  a certain  limited  range 
(see  also  Section  4.3). 
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5.6.2  Effect  of  Component  Cost  Var iat  ion s 

We  examine  the  fo.ll owing  three  cost  components:  earth  station 

cost,  space  segment  cost  and  terrestrial,  mileage  cost.  Figure  10 
shows  that  the  effect  of  unit  component  cost  changes  on  the  optimal 
satellite  switch  number.  We  observe  that  changes  of  either  ter- 
restrial mileage  cost  or  earth  station  cost  has  a more  pronounced 
impact  on  optimal  number  of  satellite  switches  than  changes  in  the 
unit  space  segment  cost.  Similar  to  results  obtained  in  [NAC,  1976] 
we  would  expect  that  the  effect  of  a unit  component  cost  variation 
on_  the  optima m_  number  of  satellite  switch e s is  rough  1 v_  propo rt ienal 
to  the  rate  of  change  of  that  component  with  respect  to  the  number 
of  satellite  switches  at  the  optimum  point.  A comparison  between 
Figures  9 and  10  also  supports  this  claim. 

5.6.3  Effect  of  Traffic  Level 

In  Figure  11,  we  plot  the  curve  of  optimum  number  of  satellite 
switches  as  a function  of  the  traffic  level,  for  various  values  of 
unit  earth  station  cost,  for  the  11  ADD-algorithm-generated  back- 
bone switch  set. 

As  expected,  for  a given  unit  earth  station  cost,  the  optimal 
number  of  satellite  switches  increases  as  the  traffic  requirement 
level  increases.  Moreover,  the  satellite  switch  selection  problem 
is  only  meaningful  over  a limited  traffic  range.  When  the  traffic 
level  is  below  this  range,  satellite  communication  is  not  cost- 
effective  (i.c.,  no  earth  stations  need  be  used),  and  when  traffic 
level  is  above  this  range,  earth  stations  should  be  used  at  all 
backbone  switch  sites. 

It  can  also  be  seen  from  Figure  11  that,  for  a given  traffic 
requirement  level,  the  number  of  satellite  switches  decreases  as 
the  unit  earth  station  cost  increases. 
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The  optimal  satellite  switch  number  versus  traffic  requirement 
curves  in  Figure  11  exhibit  an  interesting  property:  For  each 

specific  earth  station  unit  cost,  there  seems  to  be  a threshold 
traffic  level  below  which  the  number  of  feasible  satellite  switches 
increases  slowly  as  the  traffic  level  grows,  but  beyond  which  the 
number  of  feasib1  satellite  switches  increases  rapidly,  until  all 
switches  are  used  as  satellite  switches.  This  can  be  explained 
intuitively  as  follows.  For  each  system,  there  are  usually  some 
natural  traffic  concentration  clusters.  For  example,  in  the 
AUTODIN  II  system,  there  is  a cluster  in  the  West  coast  area,  a 
cluster  in  the  Central  area,  a cluster  in  the  Northeastern  area, 
and  a cluster  in  the  Southeastern  area.  Nodes  in  the  same  cluster 
are  much  closer  than  nodes  in  different  clusters.  Thus,  as 
traffic  level  increases,  one  node  from  each  cluster  quickly 
gets  selected  as  the  satellite  switch  site.  Then  for  a large 
range  of  traffic  values,  very  few  additional  nodes  can  be 
profitably  used  as  satellite  switch  locations.  Finally  when 
the  traffic  level  passes  a certain  threshold  value,  even  nodes 
in  the  same  cluster  can  communicate  to  each  other  effectively 
through  satellite,  thus  many  nodes  become  attractive  for  earth 
station  sites  almost  simultaneously. 

5.7  Procedure  Run  Time 


The  run  time  for  the  MODULARIZED  ADD  algorithm  has  been  re- 
ported in  [NAC,  1976].  For  the  present  application,  since  we  are 
not  doing  detailed  backbone  cost  tradeoff  evaluations  (see  Section 
3),  the  run  time  should  actually  be  several  times  faster. 

The  run  times  for  the  satellite  switch  selection  algorithm 
are  recorded  in  Table  5.  As  expected,  the  algorithm  is  quite 
slow.  Thus,  for  switch  sets  with  more  than  15  switches,  some  of 
the  satellite  switches  have  to  be  preselected  and  some  of  the 
switches  have  to  be  excluded  in  order  to  have  a reasonable  run 
time. 
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i 


1 

SWITCH  SET  SIZE 

RUN  TIME* 

1 ^ 

8 

O'  5" 

1 

11 

O'  22" 

, 

16 

15' 44” 

TABLE  5:  CPU  TIMES  VS.  NUMBER  OF  SWITCHES 

SATELLITE  SWITCH  SELECTION  ALGORITHM 

Computer  Used:  PDP-10 

*0'5"  means  0 minutes  and  5 seconds. 
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6.  CONCLUSIONS 

In  this  chapter,  thr  problem  of  optimally  locating  the  terres- 
trial and  satellite  switches  in  a large  scale  network  environment 
is  concisorcd.  This  problem  is  generally  recognized  to  be  an  im- 
portant factor  in  the  overall  optim.*  n ol  packet  switched  net- 

work designs  [GERLA,  1974b] . 

Two  procedures  are  presented  The  first  procedure  selects 
the  backbone  switches  from  a set  of  candidate  switches  for  a mixed 
media  network.  This  prr - adure  is  a simple  modification  of  the 
MODULARIZED  ADD  algorithm  presented  in  [NAC,  1976]  and  is  very 
efficient.  The  second  procedure  selects  the  satellite  switches 
from  the  backbone  switches.  This  procedure  is  basically  ex- 
ponential . 

Mixed  network  designs  based  on  the  switch  locations  selected 
by  the  above  two  procedures  were  generated.  It  is  found  that  mixed 
network  designs  generally  offer  considerable  savings  over  network 
designs  using  only  terrestrial  technologies.  Moreover,  it  was  ob- 
served that  the  number  of  satellite  switches  selected  is  generally 
quite  • ^dependent  of  the  backbone  switch  set  size.  Hence,  for  a 
large  switch  set,  the  selection  can  be  restricted  to  switch  sets 
of  size  within  a certain  limited  range. 

Parametric  studies  were  carried  out  to  determine  the  sensitiv- 
ity of  the  solutions  to  changes  in  some  of  the  network  parameters. 
It  was  found  that,  for  a given  backbone  switch  set,  the  optimum 
number  of  satellite  switches  is  more  sensitive  to  the  variations 
in  the  unit  earth  station  cost  and  terrestrial  mileage  cost,  and 
loss  sensitive  to  variations  in  the  unit  satellite  space  segment 
cost.  This  is  attributed  to  the  fact  that  the  rate  of  change  of 
the  eui Lh  station  cost  and  the  terrestrial  mileage  cost  as  a func- 
tion of  the  number  of  satellite  switches  is  higher  than  the  rate 
of  change  of  the  space  segment  cost. 


6.1 


NETWORK  ANALYSIS  CORPORATION 


Some  possible  areas  for  further  investigation  are  the  follow- 
ing : 


1.  Development  of  more  efficient  algorithms  for  satel- 
lite switch  selection  (possible  using  clustering  tech- 
niques) . This  is  especially  important  for  a large  back- 
bone switch  set. 

2.  Extension  of  the  algorithms  to  handle  more  general 
switch  models  (e.g.,  monlinear  switch  cost). 

3.  Development  of  lower  bounds  to  assess  the  accuracy 
of  heuristics. 
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APPEND] X A 

GENERATION  OF  REFLECTED  BINARY  CYCLIC  CODES  (GRAY  CODES) 

In  this  Appendix  we  derive  a procedure  for  generating  the  suc- 
cessive reflected  binary  cyclic  codes  (Gray  codes) . This  procedure 
is  essentially  the  same  as  the  one  given  in  (BCOTHROYD,  1964], 

The  reflected  binary  code  may  be  considered  as  an  unusual  nota- 
tion for  binary  counts.  Its  structure  can  be  seen  in  Table  A.l. 

Counting  in  the  reflected  binary  code,  we  start  in  the  same 
manner  as  in  the  true  binary  number  system:  0,1.  Then  forced  to 

introduce  a second  position,  we  reverse  the  count  in  the  first 
position:  11,  10.  Installing  the  third  position,  we  repeat  all 

previous  counts  in  the  first  two  positions  in  inverted  sequence: 

110,  111,  101,  100.  Following  the  same  principle,  we  obtain  the 
first  four-digit  number  1100  and  count  up  to  1000,  the  highest 
four-digit  number.  Prints  at  which  we  install  a new  position 
and  reverse  the  count  are  marked  by  a horizontal  line  in  Table  A.l. 
It  can  be  seen  that  the  reflected  binary  code  is  cyclic,  i.e.,  when 
counting  from  one  number  to  the  next,  only  one  digit  is  changed. 

Let  (9n*  ^n-l/',,#  ^1^  Gray  code  word  (where  g^ 

is  the  most  significant  digit  and  g,  is  the  least  significant 

t 

digit)  and  (bn,  ...,  b-^)  is  the  k binary  code  word.  It  can  be 
shown  that  [GSCI1WIND,  1975]; 

g.  = b-  0 b.+1,  for  i = 1,  ...,  n-1 

(A.l) 

2 addition. 

...,  g^)  and  (g  , g 1#  ...,  g^)  be  two  con- 

and  (b  , b^  ,,  ...,  b. ) and  (b  , b ,,...,  b. ) 
n n- 1 I n n- l I 


g = b , 
^n  n 


where  © denotes  mod 

Lct  (V  gn-l' 

secutive  Gray  codes 
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REFLECTED 
BINARY  CODE 

94g3g2gl 

BINARY 

EQUIVALENT 

b4b3b2bl 

0 0 0 0 

0 0 0 0 

0 0 0 1 

0 0 0 1 

0 0 11 

0 0 10 

0 0 10 

0 0 11 

0 110 

0 10  0 

0 111 

0 10  1 

0 10  1 

0 110 

0 10  0 

0 111 

i ro"o 

10  0 0 

etc . 

etc . 

TABLE  A.l:  REFLECTED  VERSUS  TRUE  BINARY  NUMBERS 
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bo  the  two  corresponding  consecutive  binary  codes.  Let  i be  the 
smallest  k such  that  = 0 (i.e.,  b^  = 0 and  b.  = ...  = bj  = 1). 
Then  it  can  be  easi  ly  shown  that 


1.  bk  = b^ , for  k ~ i + 1 , 


n 


2.  b.  = 1 

l 

and 

3.  b^  = 0 , for  k = 1,  . ..,  i-1. 


Combining  the  above  with  Formula  (A.l),  we  obtain  the  following 
relationship  between  (g)  and  (g) : 


1.  gk  = gk,  for  k - i + 1, 


n; 


2.  g.  - b.  © b.+1  = b.  © 1 © b.+1  = gA  © 1; 


and  for  i > 2 , 


3‘  gi-l  = gi-l  = 1; 

4.  gk  = gk  = °,  for  k = 1 , . . . , i - 2. ’ 


Consequently,  we  have  the  following  two  cases: 


CASE  1 i = 1. 

In  this  case,  the  pattern  is; 


bn  . . . b2o  * — gn  • • • ^2  yi 


b ...  b,l 
n 2 


>n  •••  g2  gl 


It  follows  that  g^  = b 2 and  g^  = b2  © 1 = g^  © 1- 
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CASK  2 i > 2. 


In  this  case,  the  pattern  is; 

b . . . b . ^ . 0 1 . . . 1 < ► g ...g.,g.  10...0 

n i + l Jn  ^ H 1 J i 

b ...b.,.10...0  < — > g ...  g . , . g. 

n l+l  Jn  Ji+l^i 

It  follows  that  g.  = b.,,  and  g.  » b. , , © 1 = g.  © 1, 

l+l  ^i  l+l  Ji 


i + l t ■*’ 

o 

o 

® 1 = g. 

© 1. 

Notice  that  Case  1 occurs  for  every  other  code  in  the  sequence 
(every  even  number  has  b-^  = 0).  Also,  in  Case  2,  i-1  corresponds 
to  the  smallest  k such  that  g,  = 1.  These  observations  lead  to  the 
following  Gray  code  generation  procedure.  Th. s procedure  is  essen- 
tially the  same  as  the  one  given  in  [BOOTI1ROYD,  1964  J . 


Procedure  for  Gray  Code  Generation 

1.  Initialization: 

g^  «-  0 for  i = 1,  ...,  n+1. 
E <-  1. 


2.  If  E = 1,  then  ic  l ; 

otherwise,  let  i - 1 be  the  smallest  k such  that  g^  = 1. 

3 . g . «-  g . © 1 
E «-  E ® 1. 
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4.  If  <j  j - 1,  then  Done ; 

otherwise  go  to  2. 


Notice  that  in  the  above  procedure,  the  stopping  condition  is 
effected  when  gn+^  = 1 (which  signifies  that  all  n-bit  cod^s  have 
been  enumerated) . 

Next,  wc  consider  the  time  complexity  of  the  above  procedure. 

Notice  that  for  the  n-bit  Gray  code  sequences,  there  are  2n  ^ codes 

with  E = 1,  2n  2 codes  with  g,  - 1,  E = 0,  2n  2 codes  with  g = 1, 

n — K “ 

g^  = 0,  E = 0,  and  in  general  2 codes  with  g 2 = •••  = 

g^  = 0,  E - 0.  Now,  for  a code  word  with  g^  = 1,  g.  2 • • • = 

g^  = 0,  E = 0,  the  amount  of  work  involved  in  getting  the  next  code 
is:  k-1  steps  to  find  that  9^-1  " ^ an<3  ^y-2  = ***  = y^  ” 0,  1 step 

in  setting  g^  to  g^  © 1,  and  1 step  in  changing  E to  E 0 1.  lienee, 

n 

the  total  amount  of  computation  involved  is  1.  (i  + 1)  2n  1. 

i-1 

n - * +1 

In  the  following  assertion,  we  show  that  }]  i 2n  1 = 2n  -n-2. 

i.=o 

Hence,  the  total  computation  is  essentially  2n+^  + 2n. 

ASSERTION 


^ i ~n+l  ~ 

l 1 2 =2  - n-2. 

i-o 


PROOF 


1;  i 2n"1  = r.  n2n_1  - r.  (n-i)  2n_1 


n 

X 

i=o 


1=0 


1=0 


n n 

= n 1 2 X i2 

i~o  i=o 


i~o 


n 


i=o 


(I) 
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n . n+1 

..  1 x 

i\ow  / >.  x = £ 

i-o 

Thus,  we  have 


i;  2 1 = 2n+1  - l. 

i=o 


Moreover,  differentiating  (II),  we  obtain 


v • i~l  n+1  / , , \ n_._i 

l ix  = nx  - (n+  J ) x +1 . 


i=l 


(x-1) 


(ID 


(HI) 


Consequently , 

>;  i21-1  = n2n+1  - (n+1)  2n  + 1, 

i=l 

i.e.  , 

1 121  = n2n+2  - (n+1)  2n+1  + 2.  (IV) 

i-o 

Combining  (I) , (III)  and  (IV) , we  obtain 

l i2n_1  = n(2n+1  - 1)  - [n  2n+2  - (n+1)  2n+i  + 2) 
i=o 

= 2n+1  - n-2.  Q.E.D 


Attached  is  a copy  of  the  Gray  code  generation  procedure 
written  in  RATFOR,  and  a sample  output  with  n = 5. 
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APPENDIX  _H 

A BOUND  ON  Till;  SATELLITE  SWITCH  SELECTION 


Here  we  discuss  several  other  possibilities  in  reducing  the 
problem  size  of  the  satellite  switch  selection  problem  presented 
in  Section  4. 

Given  two  switch  locations,  A and  A'.  Suppose  that  a satel- 
lite switch  has  been  implemented  at  location  A'.  Intuitively,  if 
location  A is  too  close  to  A’,  it  will  not  be  cost-effective  to 
implement  a satellite  switch  at  A.  In  the  following,  we  derive 
a bound  on  this  "closeness". 

Take  any  switch  location  B.  Let  S(B)  be  the  satellite  switch 
for  B (S (B)  may  be  3 itself) . 

We  consider  the  possible  (terrestrial  and/or  satellite)  line 
saving  of  adding  a satellite  switch  at  A. 

Let 


C,  t-,  = Cost  of  transmitting  a unit  of  flow  from  A to  B , 

A,B  J 

assuming  no  satellite  switch  at  A, 


S (h)  = Nearest  satellite  switch  to  A, 


C^  ^ = Cost  of  transmitting  a unit  of  flow  from  A to  B , 
assuming  a satellite  switch  at  A. 


Then , 


„TER  rTER 

^A,S(A)  - ^A,A' ' 


_ . ,„TER  , ^TER  , „ 'TER, 

CA,B  mln  {CA,S(A)  CB,S(B)  CS'  CA,B}' 
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and 


'A , B 


= mm 


{ TKR 
lUB,S(B) 


+ CV 


A,  13 


If 


PTER  . rTER 
CA,S(A)  H3,S(B) 


+ 


,TER 
A , B ' 


then 


_ TER  CTER  c 

CA ,fi  ~ CAfS(A)  B,S(B)  S 


< c?ET  + c™- 

~ A, A*  B, 


+ cs' 


and 


C = CTER  4-  C . 
CA,B  CB,S(B)  + ^S’ 


Hence,  it  follows  that 


- < ^TER 

CA,B  ~ CA,B  - ^A,A’  • 


On  the  other  hand,  if 


TER  CTER  + c > cTER 

CA,S(A)  CB,S(B)  S A,B 


then 


C 


A,  B 


,TER 

A,B* 
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In  this  case,  il 


t = rTKR 
A ,13  A # H # 


then 


'A,B  “ CA,B 


rTER 
- A, A' 


o < c; 


and  if 


C 


A , B 


^.TLR  £ 

B,S(B)  S' 


then 


C 


A,  B 


C = CTEH  - (CTER  +C  ) 
,T,B  A,B  ^ B,S(B) 


< (CTER  +CTER  +C  ) 
* A,S(A)  ^B,S(B)  ^S} 


(CTER  +C  ) 
‘ B,S(B)+<V 


=.  c 


TER 

A,S(A) 


< CTER 
- CA,A’ * 


We  thus  conclude  that  in  all  cases, 


C - C < cTER 
A,B  A, B - A, A' * 


It  follows  that  the  total  saving  of  implementing  a satellite 
switch  at  A,  SAVE,  is  bounded  by; 
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J 


SAVE  = (lino  saving)  - (ground  station  cost) 


^ tA,D(CA,B-i'A,Ii)  ' _ CGRD 


i Vb’O  - CGRD 


TER 

= t v c 1 - r 

A A, A*  ^GRD' 


where 


= Total  lov;  priority  traffic  requirement  at  switch  A , 

tA  0 = Low  priority  traffic  requirement  between  A and  B. 

We  tnus  have  trie  followincj  result: 

ASSERTION 

Give  two  switches  A and  A’,  suppose  A’  has  a satellite  switch. 
If 


t 


A 


x C 


TER 
A,  A* 


< r 

- GRD 


(B.  1 ) 


(where  t^  is  the  total  low  priority  traffic  requirement  at  A) , then 
it  is  not  cost-effective  to  implement  a satellite  switch  at  A. 

Other  problem  size  reduction  bounds  are  possible.  For  example, 
suppose  a switch  is  of  some  distance  away  from  all  other  switches. 
Then  most  likely  it  would  be  cost-effective  to  implement  a satellite 
switch  at  that  location.  However,  the  bound  in  this  case  appears  to 
be  not  as  good  as  that  given  in  Eq. (B.l).  Reduction  consideration 
for  a cluster  of  switches  (e.g.,  at  at  most  two  satellite  switches 
in  a group  of  four  switches)  can  also  be  developed  along  the  same 
line  as  that  for  a single  switch. 
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CHAPTER  2 

AN  ALGOR  I TliM  FOR  DESIGN  OF  NON’ II  1 ERARCII T CA1.  C I ECU  I T-SW  ! TC  HE!)  NETWORKS 


1 • INTRODUCTION 

Over  the  last  few  months,  attention  has  been  directed  at  de- 
veloping and  proqramming  an  algorithm  tor  accomplishing  the  design 
of  a circuit  switched  network.  In  particular,  effort  has  been 
focused  at  developing  a program  for  designing  a circuit  switched 
network  which  is  non-hierarchical  and  has  a topology  based  upon  some 
"minimum  cost"  criterion.  It  is  desired  that  the  network  produced 
be  capable  of  adequately  handling  integrated  (mixed  voice  and  data) 
traffic.  Circuit  switched  network  design  for  integrated  traffic  is, 
of  course,  of  interest  by  itself.  However,  in  addition,  when  the 
wo: k in  this  area  is  coupled  with  our  work  on  packet  switched  network 
design  for  mixed  voice  and  data,  insight  can  be  gained  for  consider- 
ing the  network  design  problem  when  integrated  traffic  is  to  be 
handled  by  hybrid  switching  techniques. 

This  chapter  reports  on  the  development,  to  date,  of  a circuit 
switched  network  design  program.  The  work  carried  out  is  summarized 
in  the  sequel  as  follows.  Section  2 describes  the  "ideal"  structure 
of  a program  fer  carrying  out  circuit  switched  network  design. 

Section  3 deals  with  the  actual  programmed  algorithm  which  has  been 
developed.  It  describes  the  structure  of  the  design  program  and 
discusses  the  philosophy  of  the  design  approach.  Experiments  that 
can  be  performed  with  the  programmed  algorithm  in  the  near  term  are 
suggested  in  Section  4.  The  future  development  of  the  program  is 
discussed  in  Section  5. 
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2.  "IDKAI,"  PROGRAM  GTRUCTUR1! 

IL  is  worthwhile  to  describe  the  "ideal"  structure  of  a pro- 
gram for  designing  a circuit  switched  network.  Such  a description 
serves  two  purposes: 

1.  It  provides  a standard  against  which  to  judge 
the  structure  of  the  actual  programmed  algorithm 
which  lias  been  developed, 

2.  It  provides  a model  for  use  in  developing  design 
algorithms  in  the  future. 

Any  procedure  for  carrying  out  a network  design  must  begin 
with  certain  inputs.  In  the  context  of  the  circuit  switched  net- 
work design  under  consideration  these  inputs  are  taken  to  be: 

1.  The  set  of  nodes,  including  their  locations, 

2.  The  traffic  requirements;  the  node-to-node  average 
traffic  to  be  handled  by  the  network, 

3.  A set  of  performance  requirements  to  be  satisfied 
by  the  design;  this  will  almost  always  include  some 
measure  of  "loss  probability." 

Ideally,  the  goal  of  a network  design  procedure  is  to  specify  a 
minimum  cost  network  of  the  input  node  set  which  satisfies  the 
traffic  and  performance  requirements.  Specification  of  a network 
of  nodes  is  taken  to  mean: 
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1.  Refining  the  connectivity;  listing  the  node  pairs 
which  arc  connected  by  (directed  or  undirected)  links, 

2.  Specifying  the  routing  for  the  connectivity;  listing 
the  sequence  of  links  traversed  by  traffic  from  a given 
source  node  to  a given  destination  node, 

3.  Sizing  the  links  - for  each  link  listing  the  number 

of  trunks  required  in  the  bundle  - "the  trunk  engineering." 

Cost,  traffic  and  performance  are  related  in  a very  complicated 
and  extremely  non-linear  manner.  Except  in  the  simplest  cases,  it 
is  all  but  impossible  tq  obtain  an  optimum  procedure  for  designing  a 
minimum  cost  network  . Instead,  design  procedures  are  heuristic 
with  the  goal  of  producing  a network  which  satisfies  traffic  and 
performance  requirements  in  an  "economic"  fashion  - a network  which 
is  close  to  minimum  in  cost.  Design  procedures  tend  also  to  be 
iterative.  Given  inputs,  an  initial  pass  at  design  is  made.  This 
first  pass  design  is  then  sequentially  altered  until  the  "designer" 
is  satisfied  with  the  resulting  network  cost  and  is  convinced  that 
traffic  and  performance  requirements  have  been  met. 

The  "ideal"  structure  of  a program  implementation  of  a general 
"heuristic"  circuit  switched  network  design  algorithm  should  be 
similar  to  that  illustrated  in  Figure  1.  Here,  "connectivity,'' 
"routing,"  "trunk  engineering,"  "costing,"  and  "performance  analysis," 
are  kept  as  separate  and  distinct  modules  in  the  program.  Thus,  the 
functional  operation  of  each  of  these  elements  of  the  design  procedure 
can  be  altered  without  affecting  the  other  elements.  This  allows 
experimentation  with  the  overall  design  procedure  to  be  performed 
easily.  For  unple,  by  keeping  all  modules  the  same,  but  by  chang- 
ing the  function  code  in  the  costing  module,  the  variation  of  final 
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network  desiqn  with  different  tariff  strucLircs  can  be  observed.  By 
keepi ng  all  modules  the  same,  but  changing  the  routing  scheme,  the 
effect  on  overall  performance  can  be  judged.  The  boxes  labeled 
"instructed  adjustments,"  are  the  variations  ordered  in  the  con- 
nectivity, routing  and  trunk  engineering  modules  in  order  to  allow 
cost  to  be  reduced  and/or  performance  requirements  to  be  met.  These 
adjustments  might  either  be  automated  as  program  code  (modules 
themselves)  or  might  represent  the  effect  of  a "human  designer" 
interacting  with  the  programmed  design  procedure. 
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3.  TUP,  CURRENT  CIRCUIT  SW  l TCHKP  NETWORK  PR()OI^\M 

The  circuit  switched  network  design  [urogram  which  has  been 
developed  is  eased  upon  an  algorithm  of  J.H.  Knepley  ( KN KPLl'Y , 

1973).  There  are  other  approaches  to  the  problem  of  circuit  switched 
network  design,  e.g.,  [COVO,  1973),  [KATZ,  1973).  However,  Knepley' s 
procedure  provides  an  attractive  starting  point  for  a number  of 
reasons.  It  is  very  general,  dealing  with  a wide  variety  of  design 
issues  ranging  from  routing  to  trunk  engineering.  It  has  the  flexi- 
bility to  deal  with  many  interesting  design  parameters  such  as 
tariff  structures  and  switch  cost.  Finally,  the  design  problem  that 
the  procedure  addresses  is  that  dealt  with  in  the  consideration  of 
such  existing  networks  as  AT&T  and  AUTOVON . 

A block  diagram  of  the  developed  network  design  program  is 
illustrated  in  Figure  2.  Henceforth,  this  t rogram  will  bo  referred 
to  as  "the  baseline  program."  Although  oui  ultimate  goal  is  tne 
handling  of  integrated  traffic,  to  date  the  function  code  of  indi- 
vidual blocks  in  the  program  has  been  generated  considering  only 
voice  (or  sta t i st i cal ly  equivalent)  voice  users.  In  the  figure,  the 
correspondence  has  been  noted  between  blocks  (or  groups  of  blocks) 
of  the  baseline  algorithm  and  blocks  in  the  ideal  structure  of 
Figure  1.  Observing  Figure  2,  the  baseline  algorithm  consists, 
mainly,  of  four  segments.  The  first  is  an  "initialization  segment." 
Here  variables  are  defined,  data  read  in  and  parameters  arj  initial- 
ized. It  ends  with  the  "initialization  box."  The  second  is  a "design 
segment."  This  extends  fiom  the  box  labeled  iteration  control  to 
the  "cost  decision  box."  This  segment  is  the  "heart"  of  the  desi  n 
procedure.  Connectivity,  primary  routing  and  trunk  sizing  are 
performed  here.  The  third  segment,  is  the  alee  mate  routing  module. 

At  this  point  an  alternate  routing  plan  is  specified  for  "the  partial 
design"  which  has  just  been  obtained.  The  last  segment  is  the 
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network  performance  analysis  module.  In  this  element  the  performance 
of  the  complete  network  design  (connectivity,  primary  and  alternate 
routing  plans,  trunk  sizing  ) is  determined,  the  averaqe  node-to- 
nodc  loss  probability  is  estimated.  At  present  there  is  no  automatic 
comparison  in  order  to  judge  whether  perfoi  manco  specifications  have 
been  met  or  have  been  met,  with  too  much  maiqin.  There  is  no  feedback 
from  this  module  to  the  design  procedure. 

A description  of  the  operations  carried  out  in  each  of  the  four 
segments  of  the  program  will  now  be  given. 

3 . 1 I n i t i a lization  Segment 

In  this  segmeni  required  constants  (such  as  N,  the  number  of 
nodes)  are  stored.  Input  data  consisting  of  1)  the  node  labels  with 
the  latitude  and  longitude  of  each  node  and  2)  the  node- to-nocle 
traffic  requirements  are  supplied  in  matrix  form.  From  this  input 
data  the  internode  distances  are  computed  and  stored.  To  date, 
although  this  will  change,  the  internodc  distance  has  been  considered 
to  be  "true"  link  cost  per  trunk  between  two  nodes.  The  internode 
distance  matrix  is  equivalently  a "true"  cost  matrix.  Henceforth, 
it  is  called  C (I , J) . 

At  the  end  of  this  segment  parameters  are  initialized  for  use 

in  the  subsequent  "design  segment."  The  "current"  total  cost  of  the 

20 

network  is  set  initially  at  infinity  (10  in  the  program).  The 
procedure  in  the  design  segment  requires  an  a priori  value  of  link 
blocking  probability.  This  is  "E"  and  it  is  initialized  here  at 
some  appropriate  value,  for  example,  at  0.005.  In  later  program 
development  it  may  be  more  convenient  to  read  E from  an  input  data 
file.  The  design  procedure  finds  minimum  cost  routes  based  upon  a 
marginal  link  cost  per  trunk  rather  than  the  true  link  co^t  per 
trunk.  The  marginal  link  cost  per  trunk,  for  a link  between  I and  J 
is  obtained  using  W(I,J),  a link  weighting.  In  the  initialization 
block  W ( I , J ) is  set  equal  to  1 for  all  I,J.  These  "weights"  are 
stored  in  a weight  matrix. 
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3 , 2 Do s i_qj  1 Segmen  t 

It  is  worthwhile  first,  to  describe,  "in  word!!,  " the  design 
philosophy  and  total  operation  of  this  segment  of  the  program.  For 
the  interested  reader  more  detailed  discussion  of  each  of  the  program 
modules  in  this  segment  will  then  follow. 

3.2.1  Overview 

The  design  segment  uses  an  iterative  procedure  in  order  to 
obtain  a "minimum  cost"  network  connectivity,  primary  routing  and 
trunk  sizing  for  the  input  node  set  and  traffic  requirements.  The 
first  iteration  begins  by  connecting  nodes  together  by  paths  com- 
posed of  directed  links.  Specifically,  each  node  pair,  ( .1  , J ) , is 
connected  by  the  set  of  directed  links  which  constitute  the  "least 
cost"  path  between  node  1 and  none  J.  Cost  in  this  first  iteration 
being  the  "true"  link  cost  per  trunk  supplied  in  the  previous  seg- 
ment. This  "least  cost"  path  is  temporarily  designated  as  the 
primary  route  between  I and  J.  Thus,  the  design  algorithm  in  this 
segment  merges  together  the  "connectivity"  and  "routing"  elements  of 
the  network  design  procedure.  Using  the  input  node-to-node  traffic 
requirements  and  the  temporarily  designated  primary  routing,  each 
link  of  the  "current  network"  (the  network  of  this  iteration)  is 
loaded  with  traffic.  The  Erlang  B formula  in  (1)  is  then  used  to 
determine  the  number  of  trunks  required  by  each  link  in  order  to 
satisfy  the  a priori  designated  link  blocking  probability,  E with 
this  loaded  traffic. 

5! 

o = N! 

1 + G + G2  + G3  + . . . + GN 


2! 


3! 
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Hero  G is  the  traffic  in  Erlangs  offered  to  a group  of  N trunks. 

"13"  is  the  probability  that  an  incoming  call  to  such  a trunk  group 
will  be  blocked.  The  total  ''true  cost"  of  the  resulting  network 
(connectivity,  routing  and  trunk  sizes  have  at  this  point  been 
specified)  is  then  computed  and  stored.  In  the  last  part  of  the 
iteration  a pew  "marginal  cObL"  is  computed  for  each  link  (i.e., 
each  node  \ lir) . This  marginal  cost  is  equal  to  the  "true  cost" 
times  a weighting  factor.  The  weighting  factor,  W^,  for  node  pair 
(i,j)  is  given  by 


8E 

s , a / 3 a 
3E 

s , a/3  s 


“s , a 


= 13, 


S = 3 . . 


a=a . . , i/i 
i j / J 


W . . =oo  i=  j 
13  J 


Here,  E is  the  Erlang  V formula,  E is  the  a priori  designated 
s , a 

value  for  the  link  blocking  probability,  s^  is  the  number  of  trunks 
computed  in  this  iteration  for  link  (i,j)  and  is  the  traffic 

loading  on  link  (i/j)  computed  in  this  iteration.  With  this  weight- 
ing the  marginal  cost  effectively  considers  links  with  large  trunk 
group  sizes  to  be  more  economical  than  links  wi th  small  trunk  group 
sizes.  The  second  iteration  is  then  begun.  The  procedure  is  the 
same  as  with  the  first  iteration  but  with  the  "least  cost"  paths 
based  upon  the  new  marginal  costs  rather  than  the  "true  costs." 
Subsequent  iterations  continue  in  the  same  manner.  At  the  end  of 
each,  a network  design  (connectivity,  primary  routing  and  trunk 
sizes)  is  produced,  the  true  cost  of  the  network  design  is  computed 
and  the  link  cost  weighting  factors  to  be  used  in  the  next  iteration 
are  determined.  The  iteration  loop  terminates  when  the  computed 
network  cost  begins  to  change  only  by  an  insubstantial  amount.  The 
final  (partial)  network  design  produced  by  the  design  segment  is 
that  corresponding  to  the  last  iteration. 
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3.2.2  Detailed  Descript i on 

A detailed  description  of  each  of  the  modules  which  implement 
the  design  segment  will  now  be  given. 

3 . 2 . 2 . 1 Iteration  Control 

This  block  keeps  track  of  the  number  of  iterations  that  the 
design  segment  has  gone  through.  It  is  entered  at  the  beginning 
from  the  initialization  block  and  subsequently  by  looping  from  the 
"Cost  decision  box. " 

3 . 2 . 2 . 2 Development  of  Marginal  Cost Matrix  Block 

Here  the  marginal  cost  matrix  for  each  particular  iteration  is 
constructed.  The  construction  is  carried  out  in  two  steps.  In  the 
first  step  the  matrix  entries  are  filled  using  the  following  formula 

MC  = C W 
"^IJ  ^ TIJ 

Cjj  is  the  "true"  cost  per  trunk  of  link  (I,J)  and  Vi1J  is  the  cur- 
rent weighting  of  link  (I,J) . C^j  and  the  first  iteration  value  of 
Wy  T are  obtained  in  the  initialization  segment.  In  the  second 

MCt  , is  set  equal  to  infinity  for  all  I=J.  This  prevents  dcvelop- 
i J 

ment  of  minimum  cost  routes  having  self  loops  at  nodes. 

3 . 2 . 2 • 3 Development  of  Primary  Routing  and  Route  Costing  Block 

In  this  block  the  "primary  routing"  is  accomplished.  For  each 
node  pair  I,J  where  I^J,  the  least  marginal  cos c path  between  I and 
J is  determined.  That  is,  the  least  cost  path  between  I and  j is 
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obtained  whore  the  link  costs  per  trunk  used  are  the  marginal  costs 
computed  in  the  previous  block.  Floyd's  algorithm  for  shortest 
paths  is  currently  used  for  determining  the  "least  marginal  cost 
pata. " This  may  be  replaced  in  the  future  by  a more  efficient 
algorithm.  This  routing  establishes  a connectivity  of  the  nodes 
(i.e.#  tells  where  links  should  be  placed)  and  specifics  the  primary 
routes  or  first  routes  which  an  incoming  call  will  attempt  to  reserve 
in  establishing  an  originating  node  - destination  node  connection. 

The  current  (for  the  iteration)  marginal  link  costs  per  trunk  corre- 
sponding to  these  primary  routes  are  also  stored. 


3 . 2 . 2 . 4 Network  Loading  Block 

At  this  point  each  of  the  links  in  the  node  connectivity  defined 
by  the  previous  block  is  considered.  Using  the  input  end-to-end 
traffic  requirements,  the  total,  traffic  offered  to  each  link  by  all 
the  primary  routes  which  utilize  the  link  is  determined  and  stored. 


3.2.2. 5 Trunk  Group  Sizing  Block 


Here  a computation  is  made  of  the  minimum  number  of  trunks 
required  by  each  of  the  links  (in  the  current  node  connectivity)  in 
order  to  satisfy  the  a priori  designated  value  of  the  link  blocking 
probability,  E.  For  a given  link  having  an  offered  load  of  A erlangs 
as  determined  by  the  previous  NETWORK  LOADING  BLOCK,  the  minimum 
number  of  trunks  required,  S,  is  determined  recursively  from  the 
following  relation 


E 


s,  A 


A • E 
s + A 


(s-1 ) , A 

e(s-1),a 
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This  relation  is  derived  from  the  Erlang  B Loss  formula.  hero, 

}•;  = ] . The  desired  minimum  value,  s,  is  the  smallest  value  of 

O ,  A 

s such  that  E * < E.  The  trunk  size  outputs  are  stored  in  a 
s,A  - 

matrix. 

3 . 2 . 2 . 6 Developmen t of  Next  Weight  Matr j x Block 

In  this  block  the  link  weighting  factors,  {W^},  to  be  used 
in  the  next  iteration  of  the  design  segment  are  determined.  These 
are  computed  from  the  linx  loading  and  trunk  sizing  of  the  "design" 
established  in  the  current  iteration.  These  new  link  weighting 
factors  allow  the  computation  of  new  marginal  link  costs  in  the 
next  iteration.  The  weighting  factor,  , corresponding  to  link 
(i,j)  is  computed  using  the  following  formula  which  is  equivalent 

to 

3Es,a/3a 

w = a(]-E)~s 

ij  a [In  (a) -E-iMs"  j 3 

Here , 


s=s . . 

ID 


a=a  . . 
ID 


1.  is  the  number  of  trunks  currently  required 
by  link  (i,j)  as  determined  by  the  previous  trunk 
sizing  block, 

2.  is  the  traffic  offered  to  link  (i,j)  as  determined 
by  the  previous  network  loading  block, 


3.  s"  is  the  integer  portion  of  s^  and  ill's") 
fined  as  follows 


is  de- 
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i|i(s" ) = -y  when  s"  1 
s"-I 

= -Y  + ^ m ^ when  s"  >_  2 

m-1 

y=  0.5772156  ...  (Euler's  constant; 

The  outputs  cf  this  block  are  stored  in  the  weight  matrix,  W(I,J), 
replacing  the  current  entries  ol'  this  matrix.  Those  weights  make 
the  marginal  cost  smaller  for  efficient  links  utilizing  many  trunks 
and  larger  for  inefficient  links  where  few  trunks  are  employed.  Ef- 
fectively, link  replacement  will  ultimately  be  carried  out  by  this 
change  of  the  weight  matrix.  It  is  automated  in  this  design  pro- 
gram whereas  in  other  design  procedures  it  might  be  the  task  of  a 
human  designer  interacting  with  the  program. 

3. 2, 2. 7 Network  Costing  Block 

Here  the  cost  of  the  network  design  obtained  in  the  current 
iteration  is  computed.  The  following  formula  is  used  to  compute 

cost. 

► 

Cost  = C.  . T.  . 

ID  ID 
i t j 
i^j 


is  the  number  of  trunks  employed  by  link  (i,j)  in  the  current 
design . 
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3. 2. 2. 8 Cost  Decision  Block 


At  this  point  the  total  true  cost  of  the  network  produced  in 
the  current  iteration  (computed  in  the  previous  block)  is  compared 
to  the  true  cost  of  the  network  produced  in  the  last  iteration.  If 
the  cost  of  the  current  network  is  lower  (i.o.,  there  has  been  a 
cost  improvement)  then  the  program  loops  back  to  the  iteration 
control.  It  proceeds  with  the  next  iteration  using  the  weighting 
matrix  that  has  been  most  recently  computed.  If  the  cost  of  the 
current  network  is  greater  than  or  equal  to  the  cost  vjf  the  network 
in  the  previous  iteration  (i.e.,  no  cost  improvement)  the  program 
exists  from  the  design  loop.  (The  true  cost,  C,  is  initialized  at 
hence  there  will  be  at  least  two  iterations  in  the  loop.)  Upon 
exiting  from  the  loop,  the  (partial)  network  design  produced  in  the 
last  iteration  is  printed  out,  that  is,  the  matrix  of  primary  routes 
and  the  trunk  size  matrix  are  printed.  The  true  cord  of  this 
(partial)  netwoik  produced  by  the  design  segment  is  also  printed. 

3 . 3 Alternate  Routing  Module 

The  design  segment  establishes  the  primary  routing.  For  an  in- 
coming call  at  node  "i;  intended  for  node  "j",  the  network  routing 
control  first  attempts  to  reserve  the  links  on  the  primary  route 
between  i and  j in  establishing  the  end-to-end  connectivity.  If 
this  cannot  be  done  (i.e.,  if  one  of  the  links  is  blocked) , the 
call  will  be  lost  at  this  point  unless  there  is  a method  of  alternate 
routing.  This  module  provides  an  alternate  routing  plan.  For  each 
end-to-end  node  pair  it  establishes  a hierarchical  order  of  alter- 
nate routes  co  handle  overflow  traffic  from  the  primary  route. 

The  alternate  routing  module  uses  as  inputs: 

1.  The  node  connectivity  established  by  the  output  of 
the  design  segment  (the  partial  network  design  produced 
in  the  last  iteration  of  the  design  segment) , 
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2.  The  marginal  link  costs  per  1 run!;  over  tn^»  link:;  of 
the  primary  roales  of  the  dosiqn  segment  output. 


The  first  procedure  in  the  module  is  to  determine  which  of  the 
nodes  are  pendent  or  semi-pendent  relative  to  the  topolouy  established 
by  the  output  of  the  design  segment.*  The  development  of  alternate 
routes  for  traffic  from  pendent  and  semi-pendent  source  nodes  is 
handled  differently  from  the  development  of  alternate  routes  for 
traffic  from  "normal"  (all  other)  source  nodes. 

The  module  handles  the  problem  of  alternate  routing  front  the 
normal  source  nodes  first.  Consider  the  source  node  - destination 
node  pair,  i,j,  where  1 is  a normal  node.  The  program  determines 
the  set  of  nodes  {k } which  are  adjacent  to  node  i (in  the  sense 
that  there  is  a link  from  node  i to  each  of  the  nodes  in  this  set). 

The  node  in  {k } which  is  adjacent  to  i,  but  also  on  the  primary 
route  from  i to  j is  purged  from  this  node  sr*t.  For  each  of  the 
remaining  nodes  tne  parameter  d.,  . given  by 


d . . . = Z . . + 

lkj  lk  kj 


is  computed.  Here  Z.,  is  the  marginal  cost  per  trunk  of  the  link 

1 K 

from  i to  k , Zx  . is  the  sum  of  the  marginal  costs  per  trunk  of  the 

links  of  the  primary  route  from  k to  j.  Each  d..  . is  in  corres- 

1 K ^ 

pondence  with  a possible  alternate  route  from  i to  j.  Specifically, 
dik*j  corresP°nds  to  the  alternate  route  composed  of  the  link  (i,k*) 


followed  by  the  primary  route  from  k*  to  3.  The  set  of  d. 


ikj 


are 


* Pendent  nodes  have  only  one  node  adjacent  to  them  and  can  only  be 
entered  or  exited  using  the  link  from  that  node.  As  used  here, 
semi-pendent  nodes  may  have  more  than  one  node  adjacent  to  them, 
but  can  only  be  left  through  one  link. 
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ordered  from  lowest  to  highest.  The  alternate  route  corresponding 
to  the  lowest  d . is  considered  tihI  is  chocked  to  see  if  it  contains 

J K J 

any  loops.  If  it  does  not  it  is  accepted  as  the  first  in  a hierarchy 
of  alternate  i-to-j  routes.  The  procedure  continues  with  the  next 
lowest  d.,  etc.,  until  either  the  required  number  of  alternate 

IK] 

routes  is  obtained  or  until  there  are  no  more  entries  on  the  list  of 
d^j's.  If  a route  corresponding  to  a particular  d^  , is  found  to 
loop  it  is  not  considered  as  a viable  alternate  route  candidate  and 
is  discarded.  There  may  be  situations  where  for  a particular  i,j 
pair  all  alternate  route  candidates  considered  by  the  module  are 
found  to  loop.  The  module  notes  this  and  the  program  prints  this 
out  as  information  to  the  user. 

With  the  alternate  routing  from  the  normal  source  nodes  com- 
pleted the  module  handles  the  alternate  routing  from  the  pendent  and 
semi-pendent  nodes  in  the  following  manner.  Consider  the  source 
node  - destination  node  pair  i,j  where  i is  either  a pendent  or 
semi-pendent  node.  Let  k be  the  (only)  outward  adjacent  node  to  i. 
The  module  designates  the  first  alternate  route  from  i to  j as  being 
composed  of  the  link  (i,k)  followed  by  the  first  alternate  route 
from  k to  j . Other  alternate  routes  are  similarly  designated. 

3 . 4 Network  Performance  Analysis  Module 

3.4.1  Overview 

In  this  module  the  network  design  (primary  and  alternate  routing 
plans  and  trunk  size  specification)  has  its  performance  analyzed. 

This  module  has  two  outputs.  The  first  is  a matrix  of  estimates  of 
the  nodo-to-node  loss  probabilities.  The  second  is  an  estimate  of 
the  "average"  node-to-node  loss  probability.  This  is  a weighted 
average  of  the  matrix  entries  with  the  weights  being  the  fraction  of 
total  traffic  corresponding  to  each  node  pair. 
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The  module  requires  as  inputs, 

1.  The  primary  routing  plan  and  trunk  sizes  specified 

by  the  design  segment. 

2.  The  alternate  routing  plan. 

3.  The  node-to-nodo  traffic  requirements. 

4.  n,  the  number  of  primary  and  alternate  routes  for 

each  end-to-end  pair  (assumed  the  same  for  each  pair) . 

The  actual  computation  of  the  loss  probability  output  is  a 
complicated,  iterative  procedure,  and  is  best  explained  using  the 
block  diagram  illustrated  in  Figure  3.  The  "inputs  and  initializa- 
tion" box  contains  the  inputs  to  the  module  as  described  above. 

The  module  also  uses  two  matrices,  OVRFLO(I,J)  and  CAP(T,J).  Each 
has  dimensions  of  number  of  nodes  x number  of  nodes.  OVRFLOW(I,J) 
is  initially  set  equal  to  the  traffic  requirements  matrix.  CAP(I,J) 
is  initially  set  equal  to  the  trunk  size  matrix.  The  loop  in  the 
block  diagram  goes  through  n iterations.  In  the  first  iteration 
the  end-to-end  traffic  is  loaded  onto  the  primary  routes,  overflow 
traffic  is  determined  and  the  residual  capacity  of  the  network  is 
also  determined.  In  subsequent  iterations,  overflow  traffic  is 
systematically  loaded  onto  the  alternate  routes  until  the  hierarchy 
of  routes  is  exhausted  The  ''nd-to-end  loss  probabilities  are 
estimated  from  the  remaining  overflow  traffic.  The  detailed  opera- 
tion of  the  first  iteration  is  described  below;  a description  of 
subsequent  iterations  follows. 
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FIGURE  3:  BLOCK  DIAGRAM  OF  COMPONENTS  OF  NETWORK  PERFORMANCE  ANALYSIS  MODULE 
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.L 4.2  Traffic  Assignment  Block 


Each  link  in  the  network  design  is  dealt  with  separately. 
Consider  link  (i,j),  the  set  of  primary  paths  which  uses  this  link 
is  determined.  The  total  traffic  over  these  primary  paths  is  then 
computed  using  OVERrLO (1 , J)  as  the  traffic  requirements  matrix  (at 
this  point  it  equals  the  traffic  requirements  matrix  anyway) . Let 
this  total  traffic  be  denoted  by  SUM(i,j).  The  primary  paths  which 
use  (i,j)  can  be  designated  by  the  corresponding  source  node  - 
destination  node,  (I,J).  For  link  (i,j)  and  each  of  the  correspond- 
ing (I , J) ' s the  following  parameter  is  computed 


ID 


(i,  j) 


(I,  J) 


OVHRFLO(I,  J) 
" SUM  ( i , j ) 


CAP ( i , j ) 


This  is  the  amount  of  capacity  (number  of  trunks)  in  link  (i,j)  that 
this  network  analysis  procedure  considers  as  being  allocated  to  I-J 
traffic.  This  may  be  non- integral . However,  the  diaphontinc  problem 
is  ignored  at  this  point.  Finally,  in  this  blncK,  for  each  end-to- 
end  node  pair  I,J  the  parameter  X(1,J)  is  computed.  Here 


X (I , J)  - min  ID,.  ..  (I, J) 

(i,j)  U'3) 

the  minimum  being  taken  over  all  links  (i,j)  in  the  primary  route 
between  source  node  I and  destination  node  J.  X(I,J)  is  the  effective 
capacity  of  the  primary  route  between  I and  J since  the  1-J  traffic 
will  be  limited  by  that  link  in  its  primary  path  where  the  smallest 
number  of  trunks  have  been  allocated  to  it.  The  analysis  procedure 
considers  X(T,J)  trunks  on  each  link  in  the  primary  path  of  (I,J)  as 
being  available  for  the  initial  end-to-end  traffic  requirements. 


3.15 


NETWORK  ANALYSIS  CORPORATION 


3.1.3  Overflow  Determination  Block 


Jn  this  block,  for  each  end-to-end  node  pair  1,J,  the  amount 
of  traffic  which  is  turned  away  from  the  primary  route  (overflows 
it)  is  computed.  The  following  formula  is  used  to  compute  this 
overflow  traffic: 

Overflow  t; affix 
from  1-J  primary 
path 

Ilcrp'  E[x<l,J)0  + , overfi  iOW  ( I , J ) is  thc  r;rl£UKJ  B formula  evaluated 
with  (X(I,J)+li  trunks  and  an  offered  load  of  OVFRFLO  ( 1 , J ) . The 

symbol  [Vi  represents  the  greatest  integer  contained  in  Y.  The 
computed  quantity  "overflow/  traffic  from  1-J  primary  path"  then 
replaces  the  I J entry  of  the  matrix  OVERFLOW ( , ) . 

3.4.4  Residual  Capacity  Determination  Block 

In  this  block  a determination  is  made  of  the  network  capacity 
remaining  after  the  initial  end-to-end  traffic  is  loaded  onto  the 
primary  routes  by  the  method  described  above.  This  residual 
capacity  is  that  available  for  overflow  traffic  being  loaded  onto 
the  alternate  routes.  The  operation  of  this  block  is  quite  simple 
For  each  i,j,  CAP(i,j),  the  i-j  entry  of  the  capacity  matrix,  CAP, 
is  replaced  by 

CAP  ( i , i ) - ^ y ID..  . (I,J).  The  sum  is  taken  over  all  end-to-end 

(I,J)  U/3J 

pairs  which  have  primary  routes  using  link  (i,j). 


--  OVCRFLO  (1  , J ) K 


IX  (I.  ,J)-:  1-K  OVERFLOW  (1  ,J) 
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3.4.5  The  Dec  in  ion  B lock 

If  there  are  only  primary  routes  (i.e.,  no  alternate  routine]) 
the  loop  is  left  at  this  point.  Otherwise  iterations  will  continue 
as  described  later. 

3.4.6  Loss  Probability  Estimation  Block 

Here  each  entry  in  the  last  matrix  0VRFL0(1,J)  is  divided  by 
the  corresponding  entry  in  the  input  traffic  requirements  matrix. 

The  resulting  values  are  estimates  of  the  end-to-end  loss  proba- 
bility for  each  node  pair  This  is  then  printed  as  the  output 
of  the  entire  module. 

3.4.7  Subsequent  Itera tions 

The  second  iteration  uses  the  same  procedure  as  the  first 
with  the  exception  that  traffic  is  loaded  onto  the  first  alternate 
routes  for  each  end-to-end  node  pair  rather  than  the  primary  routes. 
Due  to  this,  definition  of  the  various  parameters  used  must  be 
appropriately  changed.  The  other  iterations  follow  in  like  manner, 
with  the  third  iteration  corresponding  to  the  second  alternate 
route,  etc.  When  the  entire  hierarchy  of  alternate  routes  is  ex- 
hausted the  loop  is  exited  and  the  loss  probabilities  estimated 
as  previously  described.  The  average  node-to-node  loss  probability 
is  then  computed. 
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I 

| 4.  NEAR  TERM  EXPivR ! ME  NTS 

i There  are  a variety  of  experiments  that  can  be  performed  in 

■ the  near  future  (next  one- to- two  months)  with  the  circuit  switched 

network  design  program  in  its  current  state  of  development.  These 
I experiments  can  be  divided  into  two  categories, 

| 1.  Preliminary  tests, 

2.  Circuit  switching  strategies. 

They  are  described  in  the  sequel.  In  the  following  paragraphs  E, 
as  before  is  the  "link  blocking  probability."  "P  " is  the 
(program  provided)  estimate  of  the  average  node- to-node  loss 
probabi 1 i ty . 

4 . 1 Preliminary  Tests 

4.1.1  Sample  N e t wo r k s 

In  this  experiment  the  existing  program  will  be  run  with 
several  different  traffic  matrices.  The  point  of  the  experiment 
is  tG  determine  how  global  traffic  patterns  affect  the  topology 
of  the  network  design  output.  One  example  of  the  results  of  this 
experiment  is  illustrated  by  Figure  4.  The  network  topology  illus- 
trated by  this  figure  was  obtained  when  the  design  program  was  run 
with  the  traffic  requirement  matrix  given  in  Table  1. 

4.1.2  Variation  of  Performance  With  K 

The  current  program  requires  E as  an  input.  In  one  test.E 

-2 

has  been  set  equal  to  (0.5)  10  in  hopes  of  obtaining  a network 

_2 

design  with  performance  specified  by  a PT  of  10  . Observations 
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TABLH  1:  Sample  Trof 

lie  Req . Matrix 

(Values  in 

V,  r 1 a n c[s ) 

Nod 

Nodes 

es  1 

2 

3 

4 

5 

6 

1- 

0.00000 

3.86600 

0.26900 

1.128600 

1.12900 

0.74500 

2- 

5. 9940 

0.00000 

2. 16000 

3.50600 

5.94800 

2.90000 

3- 

0. 2G200 

2.29100 

0. 00000 

0.43700 

0. 35000 

0.75900 

4- 

2.14500 

4.12600 

0. 48100 

0.00000 

1.  35300 

2.62100 

5- 

1.18600 

7.13300 

0. 60700 

1. 30600 

0.00000 

1.09000 

6- 

1.05000 

2.24200 

0. 68100 

0.96100 

0.84600 

0.00000 

7- 

2.16600 

3.65900 

0.01700 

1 . 02100 

1.62700 

2.38700 

8- 

7.69100 

7.35800 

1.22400 

7.92000 

4 . 78700 

7 . 63100 

7 

8 

1- 

1.58400 

6. 05700 

2- 

2.88900 

5. 3700 

3- 

2.02000 

1.20200 

4- 

1.69200 

3.73400 

5- 

1.07300 

2.92400 

6- 

1.47200 

3.52600 

7- 

0. 00000 

1.95900 

8- 

3.78700 

0.00000 
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to  date,  based  upon  some  initial  work,  indicate  that  this  has  not 
always  been  met.  In  this  experiment  program  structure  will  not  bo 
altered.  However,  final  network  performance  (as  represented  by  P ) 
and  network  cost  will  be  examined  as  K is  varied  through  the  follow- 
ing set,  (0.5)10  \ 10  ^ , (0.5)10  ^ . An  example  output  of  this 

experiment  is  illustrated  in  Figure  5.  The  curves  shown  illustrate 
P as  a function  of  E for  the  network  of  Figure  <1.  Curve  "a"  corrc- 
sponds  to  P being  taken  as  the  average  node-to-node  loss  probability 
over  all  node  pairs  in  the  network.  Curve  'bM  corresponds  to  the 
average  being  taken  over  all  those  node  pairs  where  the  source  node 
is  not  pendent  or  semi-pendent  (i.e.,  excluding  nodes  4 and  5.  Node 
5 is  semi-pendent  although  not  evident  in  The  figure.)  Node  8,  in 
Figure  4,  is  an  (non  pendent  node  type)  articulation  point.  Curve 
"c"  is  the  P,  obtained  by  averaging  over  those  nodes  in  the  left 
hand  partition  of  the  network  given  by  this  articulation  point. 

Curve  "d"  is  the  Pj  obtained  by  averaging  over  those  nodes  in  the 
right:  hand  partition  of  the  network. 

4 . 2 Circuit  Switching  Strategies 

4.2.1  Variation  of  P With  the  Degree  of  Alternate  Routing 

Li 


In  this  experiment  the  variation  of  P is  examined  as  the 

Jj 

number  of  alternate  routes,  provided  for  each  node  pair  communica- 
tion, is  increased. 

4.2.2  Performance  Penalty  for  Alternate  Routing 

As  Aaron  Kershenbaum  has  pointed  out  alternate  routing  may  pro- 
vide a mixed  blessing  [KERSHENBAUM,  1976] . Alternate  routed  traffic 
may  steal  capacity  from  primary  routed  traffic.  In  this  experiment 
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the  network  performance  analysis  module  of  the  existing  program 
will  be  changed  by  penalizing  alternate  routed  traffic.  With 
each  traffic  .loading  iteration  of  this  module  the  residual  capacity 
will  be  decreased  by  some  factor  so  as  to  account  for  the  inter- 
action of  alternate  and  primary  routed  traffic  in  a pessimistic 
way.  The  variation  of  P^  with  the  capacity  reduction  factor  will 
be  examined.  If  P is  very  sensitive  to  the  traffic  interaction 
a new  module  will  have  to  be  substituted  for  the  current  method 
of  performance  analysis. 

4.2.3  Effect  of  Making  Links  I3ila tera  1 

In  this  experiment  the  trunk  group  sizing  block  of  the  pro- 
gram is  alcered  so  that  the  output  of  the  design  segment  is  a net- 
work having  bilateral  (equal  capacity  in  both  directions)  links. 
Each  link  will  now  have  capacity  equal  to  the  maximum  of  what  it 
was  previously.  The  effect  of  this  on  network  (cost)  and  P will 
be  explored. 
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5 . fijtu k k n r- v f form  f.nt 

In  the  months  ahead  the  following  topics  will  be  addressed  in 
the  further  development  of  the  circuit  switched  network  design 
program : 

5.1  Reliability 

A program  module  will  be  created  which  analyzes  the  reliability 
of  the  network  design  output. 

5 . 2 Tar  i f f S true ture 

The  method  of  dealing  with  link  cost  per  trunk  will  be  changed 
so  as  to  account  for  switch  complexity. 

5 . 3 Alternate  Routing 

A wide  variety  of  other  alternate  routing  schemes  will  be 
considered  for  use  in  the  program.  These  will  include  "no  alternate 
routing,"  "progressive  routing,"  and  alternate  routing  with  total 
path  diversity.  Toward  this  end  an  exhaustive  search  and  classifies 
ticn  of  the  various  types  of  routing  strategies  that  can  be  used  in 
telecommunication  networking  is  being  performed. 

5 . 4 Mixed  Voice-Data  Traffic 

Elements  of  the  program  will  be  suitably  altered  to  account  for 
the  fact  that  traffic  may  be  the  desired  integrated  vojcg  and  data 
rather  than  pure  voice  (or  statistically  equivalent  voice)  which  has 
been  considered  so  far.  In  particular,  the  formula  for  computing 
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link  blocking  probability  will  have?  to  bo  changed.  In  making  tneso 
altorat  ions  a variety  of  cases  will  be  considered  in  order  to  test 
the  robustness  of  the  network  design.  Specifically,  we  will  consider 
the  following  traffic  situations: 


1. 

Pure  voice. 

2. 

Heavy  voice  - 

1 ight 

data . 

3. 

Equal  voice  - 

equal 

data . 

4. 

Light  voice  - 

heavy 

data . 

5. 

Pure  data. 

5 • 5 Network  Performance  Analysis  Module 

This  will  be  improved  in  order  to  obtain  sharper  estimates  of 
the  end-to-end  loss  probabilities. 

5 . 6 Incorporation  of  the  Node  Model 

The  baseline  program  requires  link  blocking  probability  as  an 
input.  This  is  related  to  link  loading  and  trunk  sizing  which  are 
currently  accounted  for  in  the  baseline  program.  However,  it  is 
also  inherently  related  to  the  properties  and  capabilities  of  the 
network  nodes  that  the  assumed  circuit  switch  models.  In  addition, 
link  cost  per  trunk  should  include  the  cost  of  termination  which  is 
related  to  switch  complexity.  To  date  these  items  have  been  neglected 
and  have  never  been  taken  into  account  in  great  detail  in  any  investi- 
gation of  circuit  switched  network  design  techniqes.  However, 
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future  dovo  1 opine  n t of  the  base  lino  algorithm  will  take  switch  com- 
plexity into  account.  This  will  bo  a major  study  effort.  A consider 
able  amount  of  study  of  switch  models  has  til  ready  boon  done.  Specif- 
ically, an  extensive  survey  and  modeling  effort  .is  currently  underway 
to  investigate  the  cost  performance  tradeoffs  associated  with  repre- 
sentative circuit  switch  architectures.  Of  necessity,  this  has 
entailed  an  assessment  of  most  existing  technology  including: 

Electromechanical  systems. 

Electronic  systems, 

Computer  controlled  systems, 

and  a variety  of  circuit  switching  techniques: 

Space  division, 

Time  division, 

Hybrid  space-time  division. 

The  short  term  goal  is  to  incorporate  the  relevant  characteristics 
of  many  switches  into  a single  generic  switch  model.  This  will,  of 
course,  subsume  all  ongoing  work  directed  at  specific  models.  The 
long  term  objectives  are  to  investigate  the  impact  of  a given  circuit 
switch  architecture  on  network  performance.  The  merits  of  specific 
architectures  can,  therefore,  bo  more  meaningfully  determined  in  the 
presence  of  network  constraints. 

Two  areas  of  particular  importance  in  which  the  circuit  switch 
can  influence  overall  performance  are  internal  call  blocking  and 
setup  delay.  Although  more  economical,  circuit  switch  architectures 
that  permit  internal  blocking  can  lead  to  intolerable  end-to-end 
blocking,  cost-benefit  tradeoffs  will  be  investigated.  The 
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additional  constraint  of  switch  setup  delay  will  bo  studied  in 
detail.  1 nterrel  a t i onshi  ps  between  routing  plans  (fixed,  alternate’ 
vs.  adaptive),  signal linq  techniques  (inband  vs.  common  channel), 
and  the  operating  speed  of  the  switch  will  be  examined  relative  to 
acceptable  values  of  cross-network  (setup)  delay.  In  addition  to 
existing  switch  hardware,  a hypothetical  new  generation  of  circuit 
switches  will  be  used  as  the  referrent  in  order  to  design  circuit 
switched  networks  which  can  satisfy  AUTOVON  requirem  nts.  Apprecia- 
ble cost  reductions  over  the  current  system  are  expected  to  emerqe . 

A rudimentary  effort  is  als  envisioned  that  will  address  the  inte- 
gration of  a second  traffic  type  (i.e.,  data)  on  existing  and  future 
circuit  switch  hardware  and  networks. 

5 . 7 An  Alternative  Circuit  Switch  Design  Approach 

Whenever  a new  design  program  is  developed  i t must  be  validated 
in  terms  of  alternative  models.  NAC  has  a very  effective  program 
to  design  packet  switched  networks.  The  program  employs  the  "flow 
deviation"  method  for  assigning  traffic  to  routes  [GERLA,  1973]. 
"Flow  Deviation"  is  an  approach  to  assigning  costs  in  a network  by 
using  convex  cost  functions  to  measure  link  saturation.  It  is  pro- 
posed to  use  the  same  program  to  design  circuit  switched  networks 
by  derivina  link  cost  functions  which  are  measures  of  blocking 
probability  [VAN  SLYKK,  19761. 
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CHAPTER  3 

ON  CONNECTIVITY  IN  MOBILE  PACKET  RADIO  NETWORKS 


1.  INTRODUCTION 


The  Packet  Radio  Network  (PRNET)  is  a store-and- forward  packet- 
switching  radio  system  [KAHN,  1975;  FRANK,  1975].  Functionally,  it 
includes  three  types  of  devices:  Terminal,  Repeater,  and  Station. 

The  repeater  provides  area  coverage  for  mobile  terminals  and  also 
acts  as  a relay  node.  The  station  provides  global  control  functions, 
gatewcty  functions  for  interfacing  with  other  networks,  and  initial- 
ization functions.  A centralized  hierarchical  ^outing  algorithm 
described  in  [GTTMAN,  1976]  and  imp’  'mented  in  the  experimental 
PRNET  is  assumed  in  this  chapter.  According  to  this  algorithm, 
packet  transportation  between  any  nodes  in  the  network  is  via  the 
station.  Fur  hermore,  the  station  initializes  and  periodically 
updates  repeater  parameters  for  routing. 

One  of  the  objectives  of  the  system  is  to  enable  communication 
between  mobile  devices.  Up  to  the  present,  investigations  have  been 
limited  to  the  study  of  Stationary  Packet  Radio  networks,  which  are 
partially  characterized  by  the  following  two  conditions: 

1.  Terminals,  repeaters  and  stations  are  stationary. 

2.  The  area  in  which  terminals  originate  is  fixed  and 

is  covered  by  appropriately  placing  the  hardware  at  some 

initial  time  t . 

o 

On  these  assumptions,  routing  algorithms  and  initialization 
algorithms  have  been  developed.  In  particular,  the  routing  algo- 
rithms use  information  on  the  network  configuration  to  derive  an 
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optimal  communica'  ion  path.  In  the  case  of  general  deployment,  in 
which  the  repeaters  and  the  station  are  not  initially  aware  of  the 
network  configuration,  some  system  capacity  must  be  spent  to  acquire 
such  information  and  initialize  the  devices  with  valid  parameters. 

The  Mobile  Packet  Radio  system  generalizes  the  stationary  sys- 
tem by  allowing  every  element  of  the  hardware  to  be  independently 
in  nr  tion  in  the  3-d imens ional  space.  Condit  ions  1 and  2 above 
must  be  dropped  and  substituted  with  appropriate  requirements,  to 
be  described  below.  There  is  a clear  distinction  between  the  two 
systems,  not  only  in  the  definition  and  analysis,  but  also  the 
problems  one  must  address.  For  example,  in  a Packet  Radio  Network 
(PRNET)  one  is  trying  to  determine  the  location  of  the  repeaters 
to  achieve  optimal  coverage  of  a certain  area.  A repeater  may  be 
permanently  placed  in  a precalculated  location.  In  the  Mobile 
Packet  Radio  Network  (MPRNET)  repeaters  are  assumed  in  motion. 

In  PRNET,  a particular  and  predetermined  zone  is  covered. 

In  MPRNET,  the  cohered  zone  varies  with  time.  Requirements 
as  to  the  coverage  of  a particular  area,  may,  however,  be  imposed. 

In  PRNET,  once  the  system  is  initialized,  it  will  remain  so; 
only  one  initialization  and  labeling  is  necessary  (except  if  there 
is  a high  failure  rate  of  repeaters) . 

In  MPRNET,  the  system  may  become  inoperative  due  to  motion  of 
the  devices.  Hence,  one  must  address  himself  to  the  questions  of 
how  long  the  system  is  operational,  or  how  long  one  must  wait  be- 
fore one  gets  coverage  at  a specified  location.  Consequently,  a 
frequent  MONITORING  OPERATION  is  required  by  the  station. 

The  dynamic  character  of  this  system  accentuates  the  n _.d  to 
acquire,  process  and  proliferate  network  configuration  information. 
Specifically,  the  connectivity  matrix  will  be  altered  within  a 
short  interval,  depending  on  th z mobility  parameters  of  the  system. 
Thus,  the  initialization  algorithms  developed  for  the  stationary 
packet  radio  must  be  activate  ' and  more  frequently,  in  addition  to 
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the  more  frequent  acquisition  of  data,  if  the  system  is  to  be  opera- 
tional. As  a consequence,  the  initialization  capacity  becomes  con- 
siderable, decreasing,  in  conjunction,  the  useful  information  capa- 
city. 

It  can  be  inferred  from  the  statements  above  that  the  avail- 
able algorithms  may  no  longer  be  optimal  in  a dymanic  environment, 
and  hence,  new  or  modified  algorithms  are  sought. 

A predictable  mobile  PRNET  (station  Knows  the  trajectories  of 
repeaters)  is  addressed  in  this  chapter.  The  efficiency  of  a pre- 
dictable system  can  be  improved  by  changing  some  protocols  in  the 
network,  as  follows: 

(i)  Connectivity  monitoring  via  ROP ' s can  be  reduced 
or  eliminated. 

(ii)  The  station  can  predict  connectivity,  and  hence 
labels,  and  send  new  labels  to  repeaters  at  the 
proper  times. 

(iii)  Alternatively,  the  station  may  send  to  a repeater 

a set  of  labels  and  associated  time  units  for  using 
each. 

It  is  noted  that  the  criterion  for  determining  labels  in  a 
mobile  network  should  not  be  based  on  "shortest  path"  only  but  take 
into  account  the  interval  of  time  during  which  a label  will  be  valid. 

A method  for  computing  connectivity  as  a function  of  time  and 
its  computational  complexity  are  presented.  Approximate  methods 
and  their  complexity  are  also  given.  A program  which  implements 
the  exact  method  was  developed  (not  presented  in  the  chapter) . 
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2 . MAT  PL  MAT  I CAL  I'  RA  ML  WORK 

In  this  section,  we  introduce  the  mathematical  framework  for 
the  MPRENT . 

Definition : [Initial  Operating  Region] 

Let  0(to)  be  an  open  '1-dimensional  region,  which  we  assume 
to  be  the  inside  of  a cube. 

Definition : [Hardware] 

Let  S = (S1,S2,  . . . ,-Sn) , R = tr1,r2,  . . . ,rm), 

T - { i ^ , i 2 / • « • » T i ■ ) be  sets  satisfying: 

a.  S n R - g , S n t = <p , t n R = 4> 

b.  SC  0(tQ) 

S C 0(tQ) 

T C 0(to) 

The  sets  S,  R,  and  i correspond  to  the  set  of  stations,  re- 

peaters and  terminals  of  a packet  radio  system,  respectively.  The 
requirement  in  a implies  that  each  device  fulfills  only  one  func- 
tion. 


2.1 


NETWORK  ANALYSIS  CORPORATION 


Def  i ni Lion : [Frame  of  Reference] 

Let  X , Yq  , Zc,  be  a cartesian  coor- 

1,  0(Lo)  1,  0(tQ)  "’l,  0(tQ) 

dinate  system,  with  axes  parallel  to  the  sides  of  O(t^)  and  with 
origin  at  S^. 

Note  that  the  frame  of  reference  can  be  arbitrarily  selected 
without  loss  of  generality.  However,  it  is  convenient  to  select 
a coordinate  system  with  the  origin  at  station  (packet  radio 
networks  under  consideration  are  assumed  to  have  at  least  one 
station) . The  set  of  free  nodes  is  implied  by  the  selection  of  the 
coordinate  system,  and  is  defined  as: 

Definition : l F r e e Nodes] 

A A 

Let  S = S - { S ^ } and  F = S U R U t.  The  elements  of  F ara 
called  free  nodes. 

Definition : [Zone  of  Coverage] 

Let  0 ( t ) be  the  zone  of  desired  coverage  at  time  t.  Possibly 

0(t)  = 0(t  ) * 0r  for  all  t. 
o o 

Definition : [Area  of  Coverage] 

Let  P be  the  x-y  plane. 

Then  P (t)  = 0(t)  H p will  be  called  the  terrestrial  area 
xy  xy 

of  desired  coverage. 

Definition : [Stationary  Geosystem] 

The  system  described  by  the  definitions  above  is  called  a 
Stationary  Geosystem,  and  will  be  denoted  as: 
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atJ  (S,  R,  T,  0 . O) 
u o 

Def  i n i t ion  : (Dynamic  Geosystem] 

Assign  each  free  node  f z F a velocity  vector  with  refer- 
ence to  as  follows: 

vf  (t)  = vf  (t) , with  | v (t)  | = v<-(t),  and  arg(vf(t))  - 0f(t) 

A Dynamic  Geosystem  is  the  system  described  by  the  previous  defini- 
tions and  the  above  vector  assignment.  Jt  will  be  denoted  by  = 

{S,  R,  i,  0(t),  V ( t ) } where  V ( t ) - (vf(t)  | f rF}. 

Equivalently,  specification  of  the  displacement  vector  of 
motion,  in  terms  of  the  x,  y,  and  z components,  could  be  given. 

Note  that  the  definition  requires  Gj  to  be  fixed  in  the  selected 
frame  of  reference,  but  no  generality  is  lost  by  this.  In  fact,  if 
in  relation  to  some  more  general  frame  of  reference,  we  had  asso- 
ciated Vg  for  all  g c S U R U t,  then  in  our  specific  frame  of  refer- 
ence through  S^,  we  would  have, 

v :=  C associated  with 

v = v'  - v'  associated  with  g r.  F 

9 9 S1  ^ 

Assumption  1: 

There  exists  a set  PR  = { , • » . , ?z ) of  transmission  powers 

for  repeaters  and  stations  with  a corresponding  set  hr  ~ { * • • • / 

y } of  transmission  ranges  (distances), 
z 
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Assumption  2 : 

There  exists  a set  = {P  } of  transmission  power  for  ter- 
minals, with  a corresponding  set  iiT  = {jj^}  of  transmission  range. 

We  assume  that  for  equal  levels  of  transmission  power,  equal  trans- 
mission ranges  are  achieved.  In  other  words,  a flat  terrain  is 
considered . 

The  zone  covered  by  transmission  by  an  element  f e F at  le\el 

U , is 


Bpt)  = (x  c Oft)  | df  x(t)  < u } 

where  d is  the  Eucledian  distance  in  3 dimensions  and  |i  i s a trans- 
mission range  from  the  sot  uR  or  n correspond ing  to  whether  f is 
a repeater  or  a terminal. 


Definition : 


V,k  Mk 

M (t)  = (a.  • (t)  ) , 1 < i,j  < n + ip . + L is  called  the 

1 / 3 

connectivity  matrix  when 


V 

i/ j 


(t) 


a 


^k 
i • 3 


(t) 


uk 

1 if  ficB.K(t.) 

0 otherwise 

1 if 

0 otherwise 


j < m + n 


j m + n 


The  top  entries  in  the  definition  of  a.  . read:  "f.  is  within 

the  range  of  device  j . " 
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Def in i t i on : [Status] 

By  status  of  the  geosystem  at  time  t,  o(t  ) , we  mean  its  con- 
nectivity matrix,  and  the  covered  zone  at  level  ji. 

Def ini tion : [Potentially  Operational  System] 

The  system  is  potentially  operational  at  time  t if  the  fol- 
lowing hold: 

1.  "Potentially  operational  W.Pv.T.  terminals:" 

0 ( t ) is  cover  ed  by  balls  of  radius  v centered 
on  r^,  and 


2.  "Potentially  operational  W.R.T.  repeaters:" 

Consider  the  graph  with  vertices  r^  and  edges 

between  r.  and  r . if  d < U,  ; then  the  graph  is 

1 1 r . , r . — k J 1 

i 1 

connected . 

Let  G ( t ) be  the  zone  for  which  the  system  is  potentially  oper- 
ational at  t.  Obviously,  we  would  like  0(t)C  G(t). 

Definition : [Actually  Operational  System] 

A system  is  actually  operational  if: 

1.  It  is  potentially  operational,  and 

2.  If  the  algorithms  and  parameters  used  by  de- 
vices for  packet  transportation  are  such  that  there 
is  a communication  path  between  any  pair  of  nodes. 

(The  communication  paths  need  not  be  optimal.) 
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This  definition  is  motivated  by  the  fact  that  even  though  a 
system  may  be  potentially  operational  (the  hardware  is  in  the  right 
position),  no  relabeling  has  yet  taken  place  - either  because  the 
connectivity  matrix  in  the  station  is  not  up  to  date  or  because 
the  station  has  not  yet  completely  proliferated  the  information 
to  the  repeaters. 

Definition : [Predictable  System] 

A system  is  predictable  if: 

o(t)  - q(o(t  ),  vf  (t),...,v,-  (t)  ) with  cj  known;  that,  is, 

° 1 1 A 

the  status  at  time  t can  be  determined  by  the  station  by  corn*' 
putation  alone. 

Definition : [Unpredictable  system] 

A system  is  unpredictable  if  either  q or  v's  are  not  known. 
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3 .  CONNECT 1 V l TY  PROBLEMS 


A number  of  pertinent  problems  can  bo  raised  at  this  point. 

The  remainder  of  this  chapter  is  dedicated  to  the  solution  of  various 
connectivity  issues  in  the  predictable  case.  The  following  ques- 
tions are  of  interest: 

1.  The  interval  of  time  during  which  repeater  r^  can 
communicate  with  the  station  S in  1-hop; 

2.  The  interval  during  which  repeater  r ^ can  communi- 
cate with  the  station  (in  any  number  of  hops); 

3.  The  interval  ot  time  during  which  the  entire  net- 
work is  connected; 

4.  The  degradation  history  of  the  network  (i.e.,  when 
at  most,  at  least,  or  exactly  k repeaters  are  unable  to 
communicate  with  the  station) ; and 

5.  The  enhancement  history  of  the  network  (i.e.f  when 
at  most,  etc.  k repeaters  return  within  the  range  of 
the  station) . 

3 . 1 Definitions 

1.  a.  A set  II  C U (51=  real  numbers)  is  said  to  be  con- 
nected if  for  all  x,yefl,  (l-a)x  + aye  II  for  all  0<a<_  1 . 

b.  { n i } ^ are  the  connected  components  of  11  if 

are  the  equivalence  classes  of  the  following  equivalence 
relation: 


x~y  iff  {(l-a)x  + ay  | 0<a£l}C  JI 
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2.  The  L spectrum  cl  ii  with  connected  components  l IL  } , *1.(1!),  is 
the  following  collection: 


* t ( II ) = (t. 


t . = inf  t } 

1 tr  FI  . 

i 


* t ( H ) = {«>} 


* c 11/4' 


if  II  = <J> 


The  points  of  *t(l!)  are  called  the  death  dates  of  II. 


3.  The  R spectrum  of  a set  II  with 
following  collection: 

t*  (II ) = { t | t . - sup  t } 

1 1 tLJl. 

t ^ t n ) = {-<»} 


components  UL  } , t A ( II ) , is  the 

if  !!/f 
if  n=cf» 


The  points  of  t*(II)  are  called  the  resurection  dates  of  H.  In  this 
section,  we  are  concerned  only  with  systems  having  a single  station 


4.  Let  R - (r^,r2#...#rj 

r.)  is  said  to  be  alive  at 
3 

with  the  station. 


. The  Joint  Life  Status  (JLS) , (r^,r?,. 

time  t if  every  repeater  can  communicate 


5.  The  JLS  is  said  to  experience  a order  failure,  k £ j,  if 

at  least  1:  repeaters  cannot  communicate  with  the  station. 
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6.  Let  r0  S,  mid  r\>r2 
The  set 


r . be 
J 


the  repi'a  I crs  ot  a geosystom. 


, k = (t  | d (r^ , ijJ  > \\}  f 0 ; i,  k j} 

represents  the  time  when  there  is  no  direct  connection  between 

and  r,  . 
k 

Figure  1 depicts  the  set  for  various  mobility  character- 

istics of  a l-repeater  geosystem. 


7.  a.  Let  1I(  be  the  time  intervals  such  that  repeater  r^ 
cannot  communicate  with  the  station  in  any  number  ot 
hups . 


b.  Let  n1,^1 , A1  be  the  time  intervals  when  at  least, 
exactly,  at  most  i repeaters,  respectively,  cannot  communi- 
cate with  the  station  in  any  number  of  hops. 

c.  Let  a1,  y1  be  the  time  instance  when  exactly, 

at  least,  at  most  i repeaters  come  back  into  range  of  the 
station  ("resurect") . 


We  shall  assume  that  the  displacement  vector  of  ivery  repeater  r^ 
is  given;  namely, 


ri:  ( x i (t)  , 


yi  (t)  , ^ (t)) 
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t 


“10 


t 


1 


“10 

kH+H+h t 


NETWORK  ANALYSIS  CORPORATION 


3 . 2 Mot i 'at ional  Example 


Consider  a 3 repeater  system.  We  wish  to  derive  * il , the  time 
repeater  r^  is  disconnected  from  the  station  S.  Such  an  event 
would  materialize  when  there  are  no  1-hop  connections,  there  are 
no  2-hop  connections,  and  no  3~hop  connections.  Each  event  (ex- 
cept the  first  one)  can  be  further  subdivided  into  subevents,  re- 
lating to  the  particular  path  structure.  A brief  analysis  shows 
that: 


"lon!~  12u  20)nt  13u"30'n(  12U  23U  30)n<  13U  32U' 20* 


2 3 

Similar  expressions  are  obtained  for  II  and  FI.  The  JLS  experi- 
ences a first  order  failure  when  either  r^  is  outside  the  range 
of  S,  or  r2  is  outside  the  range,  or  r^  is  outside  the  range. 
Thus, 


1 12  3 

n1  = 1nu  iiu  ii 


Simi larly , 


n2  = (1nn2inu(1nn3iiKj(2nn3n) , 

3 12  3 

ir  = nmifYii 


It  is  clear  that  the  quantities  at  hand  can  all  be  derived  from 
knowing  the  ~ j ( j+1)  distinct  sets.  A general  methodology 

is  now  presented. 
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4.  GENERAL  THEORY 

The  connectivity  quantities  desci ibod  above  can  be  obtained  as 
follows : 


Step  1: 


Derive  for  all  0<i<k<_j  by  solving  the  cartesian  equation 

1 K 


:jk  = {t  I (xi(t)  - xk { t ) ) 2 + (yi(t)  - y k ( t ) ) 2 
+ (zi(t)  - zk(t))2>  II2} 


Step  2 : 


Derive  hJI , l<h<j  as  follows 


in " A , . n . . ‘5h,iV5i,  i, u- 

k=l  l<i1,i2, . . . »ik£D  1 1/  -2 

i2  ^ il 

h * Vb 


.u  = . 0) 

k. 


K • i 


k-1 ' k-2 ' * 


. ,i. 


Namely, 

h 


11  = 


Of 

■.V  l<k<  j 
° k/h 


"hkU  kd  ,n[l<k?L< j ( hkU  =kLU  i )] 


o 


kj*h 
i^k,  h 
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Step 

3 : 

Construct 

II*1,  i - • h < j , as  follows: 

i.  n1  = 

V i». 

1 

2.  II2  = 

l<i<k<j  (i,i  n kn) 

3.  n3  = 

ia,k<h<j 

et  cetera. 

Step 

4 : 

Construct 

Qh,  l<h<j-l,  as  follows: 

Hh  = Ilh  - 

IIh+1 

<2 3 = iP 

Step 

Jj: 

Construct 

A*1,  l<h<j,  as  follows: 

A - y 
\=-l 


h Q1 


4 . 1 Properties 

i.  n1  on2  on3Z). . . iP 

For  example,  *11  nSl  C ^Jl , so  that  n 2 = U(1n  ^11  )C  'J 1II 


= ir 


2.  nn  represents  the  time  when  h or  more  repeaters 
are  outside  the  range  of  the  station;  IIh+1  represents 
the  time  when  h+1  or  more  repeaters  are  outside  the 
range  of  S;  hence,  indeed  = n*1  - n^+1. 
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3-  7 ' s are  disjoint  sots. 

o 

4-  11  ” time  when  at  least  zero  repeaters  are  outside 

O 

the  range,  so  II  - (-•«,«»)  . 

o 

Q - time  exactly  zero  repeaters  are  dead;  i.e., 
every  repeater  can  communicate  with  the  station;  it 
follows  that: 

n = ii  - n = i!  - ir  = c (n1) . 

where  C(A)  is  the  complement  of  set  A. 

o o 

6.  Q = A 


4 • 2 Complexity  Cons i derations 

We  define  complexity  to  be  the  number  of  computational  steps 
required  to  carry  out  a given  task.  In  this  context,  the  term 
step  does  not  refer  to  a specific  operation;  it  could  refer  to 
a multiplication  or  to  a union,  etc. 

Consider  a system  having  j repeaters  and  one  station.  The 
complexity  of  Step  1 is  j(j+l).  (Step  refers  to  the  number  of 
distinct  sets  to  be  obtained).  The  complexity  of  Step  2 is: 


Since  there  are  ( j y paths  of  k,  each  admitting  k possible  dis- 
connections. The  complexity  of  Step  3 is: 


j 

Z 

k-1 


1 


(step  refers  to  unions  and  intersections).  Clearly,  Step  2 is  the 
most  complex  of  the  three.  The  following  bounds  are  easily  ob- 
tained. 


( 2 j - 1 ) j ! < x<  e.j.j! 
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4 . 3 Approximate  Solut i ons 

In  view  of  the  high  computational  complexi ty  of  the  exact 
solutions,  approximate  solutions  are  sought.  Two  procedures 
are  proposed. 

Cumulative  A p proximo ting  Procedures  of  D e gree  d ( C A P : d ) 

In  this  approach,  a repeater  r^  is  declared  to  be  dead  if  it. 
is  unable  to  communicate  with  the  station  in  d hops  or  less.  It 
follows  that: 


11  k=l  l<i.,  ij  , iR<  j ( h,i1u  U ' ' ' " ,U  i k ^ 0 ' 


i2  * h 


1 3 ^ 1 2 ' 1 1 


1jK  71  1k-l'1k-2'  • • • (11 


Note  that  the  outer  intersection  encompasses  d-or- less -hop  connec- 
tions to  the  station. 


Sectional  Approximation  Procedure  of  Degree  d (SAP:d) 


In  this  approach,  a repeater  r^  is  termed  dead  unless  there 
is  a d-hop  connection  to  the  station.  Tnen, 


II  = , , , n 


..  ( 


l<i  * * * * # h 

h * h 

i3  / i2'il 


i u ^ U • • • • 

1 / -*■  1 x x 2 xd,u 


Ad  * Ld-v 


i 


1 
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Under  the  CAP:d,  the  computat ionaJ  complexity  is  reduced  from  ap- 
proximately (2j-i)j!  to  approximately  jc*.  The  CAP  approach  in 
superior  to  the  SAP  approach,  and  under  mild  conditions,  the  CAP : 2 

is  an  excellent  approximation  to  the  actual  answer.  With  this  pro- 

2 

cedure  of  complexity  j , up  to  30  repeater  systems  can  be  analyzed. 
Other  approximating  schemes  are  available. 

4.4  Resurection  Dates 


The  instance  when  an  inactive  repeater  comes  back  into  range 
is  useful  information,  in  particular,  in  a predictable  system. 

This  section  addresses  the  issue. 

It  is  clear  that  exactly  one  repeater  resurec4 s when  the 
state  chanqes  from  having  i dead  repeaters  to  i-1  ad.  Similarly, 
exactly  two  repeaters  resurect  when  the  state  changes,  from 
having  i dead  repeaters  to  i-2  dead  repeaters,  etc,  for  i>2.  Thus, 
it  must  follow  that: 


a1  = U (t*  (SJ1  + 1)  n *t  (U1)  ) 
i=0 


j -p 

uP  = u (t,(!i1+p)ntt(21)) 

i=0 


Then, 


6p  = 


U a 
k=.l 


P 1 

y = u a 

k=p 
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5.  EXAMPLES 
Example  1 : 


Let 


Rj:(«it  + 6.,  fi.t+y.) 


i = 1,  2 


let 


A (a1,b1,c1 ,dj ) = A (1 ) 
V+(al'b1'c1'dl>  = V(l) 


2a-jb^Cjd^  4- 

( u 2 

(-a  j djj  H-b^c  j ) f 

y/“  ,c^  ,d1 ) 

2 

, 2 

al 

+ bi 

now 


lt|  I (a.-a.)t  + (6±-6  j ) ] 2 + ((t'^-e.Jt  + (Y.-Yj)]2 


= [-»,  V_(i;j)]  yj  (V+(i;j),  oo] 


where 


v+(i;j)  = V+(a. 

* 


We  can  now  get  the  desired  connectivity  measures  by  the  pro 
cedure  of  Section  4. 
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Example  2: 

The  mobile  network  under  scrutiny  consists  of  5 repeaters 
moving  on  linear  trajectories,  within  a two  dimensional  plane  at 
low  constant  velocity.  Such  motion  is  completely  characterized  by 
specifying  the  four  parameters  of  Table  1,  where  the  units  of  time 
have  been  taken  to  be  hours,  and  the  units  of  velocity  MPH.  The 
motion  is  relat ivo  to  the  station . 

Figure  2 describes  graphically  the  trajectories  followed,  and 
Figure  3 depicts  the  status  of  the  system  at  three  sample  points. 
The  dashed  circle  represents  the  range  of  transmission  of  the  de- 
vices . 

The  solution  of  the  problem  is  furnished  by  a computer  program 
which  analyzes  such  dynamic  systems.  Topics  of  interest  arc: 

1.  The  interval  of  time  during  which  repeater  r^  can 
communicate  with  the  station  S in  l--hop; 

2.  The  interval  during  which  repeater  can  communi- 
cate with  the  station  (in  any  number  of  hops); 

3.  The  interval  of  time  during  whicn  the  entire  net- 
work is  connected; 

4.  The  degradation  history  of  the  network  (i.e.,  when 
at  most,  at  least,  or  exactly  k repeaters  are  unable  to 
communicate  with  the  station);  and 

5.  The  enhancement  history  of  the  network  (i.e.,  when 
at  most,  etc.  k repeaters  return  within  the  range  of 
the  station) . 

Table  2 shows  the  results  for  Items  1,  2,  3,  and  4b  above  for 
the  example  at  hand. 
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FIGURE  3:  STATUS  OF  THE  SYSTEM  AT  THREE  SAMPLE  POINTS 
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CHAPTER  4 

AN  APPROXIMATE  ANALYTICAL  MODEL  FOR  INITIALIZATION  OF 
SINGLE  HOP  PACKET  RADIO  NETWORKS 


1 • INTPODUCT ION 

Among  tno  objectives  of  the  packet  radio  system  are  fast  net- 
work deployment  and  the  capability  of  communication  among  mobile  de- 
vices (KAHN,  1975].  Combinatorial  analyses  {NAG,  1973]  and  simulation 
studies  [FRANK,  1975]  have  demonstrated  that  the  impact  of  the  routing 
algorithm  on  network  efficiency  is  of  prime  significance.  That  is, 
unlike  in  A jint-to-point  packet  switching  networks  in  which  a sophis- 
ticated routing  algorithm  (as  compared  to  simple  fixed  routing)  can 
result  in  an  average  increase  in  network  throughput  of  up  to  20%  to 
30%,  in  packet  radio  networks,  the  increase  in  throughput  may  be  by 
an  order  of  magnitude. 

Efficient  routing  algorithms  for  packet  radio  networks  have 
been  presented  in  [G1TMAN,  1976].  The  principal  idea  of  the 
hierarchical  routing  algorithm*  is  to  transform  the  broadcast  net- 
work to  a "point-to-point”  network  for  packet  transportation,  by 
assigning  repeaters  a code  for  routing,  called  Label.  The  problem 
is  that  the  connectivity  of  devices  is  changing,  and  therefore,  it 
is  necessary  to  develop  algorithms  for  dynamically  changing  the 
network  structure  (reconfiguration)  under  certain  conditions.  Ex- 
amples of  a changing  topology  are  when  the  network  is  mobile  (e.g., 
a Packet  Radio  System  for  a fleet  of  ships),  when  a repeater's  range 
changes  (e.g.,  from  battery  drainage)  or  when  repraters  fail.  A 
different  type  of  network  change  occurs  when  terminals  enter  or 
leave  the  areas  served  by  a given  repeater  or  station. 


* This  algorithm  is  presented  in  [GITMAN,  1976]  and  has  been 
implemented  in  the  packet  radio  experimental  system. 
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The  initialization  procedure  assumes  that  network  topology  is 
not  known  a priori  and  is  changing  with  time.  Hence,  the  initialize 
tion  procedure  involves  mapping  of  network  topology#  determining 
network  structure  (labels  for  repeaters),  and  transmitting  labels 
to  repeaters.  Initialization  and  connectivity  monitoring  algorithms 
which  are  optimum  in  some  sense,  were  presented  in  [NAC,  1975].  In 
this  chapter  we  postulate  a simple  initialization  scheme  and  perform 
worst  case  analysis  of  initialization  time  as  a function  system 
parameters . 

The  initialization  scheme  analyzed  assumes  the  following 
packet  types : 

ROP  - Repeater  On  Packet 

LLP  - Label  Packet 

ROP's  are  issued  by  uninitialized  repeaters  and  inform  the 
station  about  the  existence  of  the  repeater  and  the  link  it  records. 
LLP's  are  sent  from  the  station  to  repeaters.  A single  station  and 
single  hop  packet  radio  network  is  analyzed.  Realistic  models  are 
set  up  at  first;  however,  to  obtain  closed  form  solutions  we  must 
assume  worst  case  values  for  parameters,  obtaining  in  conjunction, 
worst  case  values  (or  upper  bounds)  for  initialization. 

The  closed  form  solution  for  the  initialization  time  is  a 
function  of  number  of  repeaters,  the  retransmission  time  intervals 
of  repeaters  and  station,  and  the  interference  (connectivity  of 
repeaters).  The  study,  by  numerical  methods,  of  the  minimum  in- 
itialization time  demonstrates: 

a.  A good  strategy  for  the  station  (in  the  absence 

of  information  traffic)  is  to  transmit  a label  every 

4 slots,  on  the  average.  This  value  was  demonstrated 

to  be  optimum  fot  large  variations  in  all  other  parameters. 

b.  For  low  connectivity  (interference)  of  repeaters, 

a repeater  should  transmit  a ROP  every  m-e  slots,  on  the 

average;  where  m is  the  number  of  repeaters. 
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2.  modi:  i, 


2.1  S ingle  Hop  Model  Without  Terminal  Traffic 

Figure  1 depicts  the  network  under  scrutiny. 

The  following  assumptions  are  made: 

1.  Single  station  with  m repeaters,  physically  at  1-hop. 

2.  No  repeater  accepts  ROP's  from  another  repeater. 

3.  Each  repeater  can  interfere  with  1 neighbors  (in 
addition  to  itself). 

4.  Common  channel,  slotted  ALOHA. 

5.  Protocol. 

a.  Each  repeater  broadcasts  ROP's,  such  that 
the  number  of  slots  between  any  two  trans- 
missions is  uniformly  distributed  on  [0,2^]. 

b.  If  a label  is  received  by  the  repeater, 
it  forwards  one  ETE  ACK  to  the  station  and  it 
halts  its  ROP  emission  except  as  in  (c)  below. 

c.  If  a repeater  has  received  a label,  it 

forwards  monitoring  ROP's  uniformly  on  [0,2(|>] 
slots,  <p  > > (Wc  shall  disregard  this 


2.1 
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FIGURE  Is  1-HOP  NETWORK 


2.2 
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effect  here,  since  by  taking  f large  we*  can 
in  effect  say  that  the  repeater  in  question 
is  essentially  quiet  after  initialization.) 

d.  If  the  station  has  issued  a label  and 
it  gets  no  ETE  ACK,  it  retransmits  the  label 
at  an  appropriate  time  - see  below. 

c.  The  station  maintains  a queue  of  all  the 
unacknowledged  labels  to  be  (re) transmitted . 

f.  Every  time  the  station  receives  a ROP 
it  checks  the  queue  of  labels.  If  there 
already  is  a label  scheduled  for  transmission 
to  the  repeater  in  question,  the  station  dis- 
regards this  ROP;  otherwise  it  adds  the  label 
to  the  queue. 

g.  FIFO  for  station  queue:  zero  processing 
and  propagation  time;  infinite  buffer  size  at 
the  station. 

2 . 2 Machinery 

Let  R^ , R2 , Rm  be  the  repeaters  in  question,  and  S the 

station . 

Let  A.  = Time  when  the  first  ROP  from  R.  reaches  S. 

1 l 

Let  Q(t)  = Those  repeaters  whose  label  is  in  the  sta- 
tion queuf*  at  time  t. 
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Q(t)  represents  those  repealers  such  that  a HOP  has  been  received 
by  the  station  by  time  L,  but  no  F.TK  ACK  has  been  received, 

A^  is  the  time  when  the  label . for  i enters  the  queue  for  the  first 
time. 

Let  B ^ = Time  when  the  label  for  i is  deleted  from  the  queue 
(namely  it  is  received  by  and  acknowledged) . 

Definition ; By  initialization  time  we  mean: 

i = max  E (B . ) ( \ ) 

i 


Let  0^  = 1 Q ( t ) | Number  of  labels  to  be  issued  at  time  t. 

We  call  0 £ the  cycle  length  at  t. 

We  now  expand  assumption  5d : 

(5d')  At  the  end  of  each  cycle,  say  t*,  S scans  Q(t*); 
the  labels  found  in  the  queue  are  then  scheduled  for 
transmission  with  an  intertransmission  delay  distributed 
U[0, 2w] . 


The  remainder  of  this  chapter  is  dedicated  to  expected  value 
computations;  hence  wo  say  that,  on  the  average,  one  ROP  every  rd 
slots  is  issued  by  a non-initialized  R^ , and  one  LLP  every  w slots 
is  broadcasted  by  S.  With  this  in  mind,  we  see  that  after  S sends 
a label  L^  to  R^,  it  times  out  for  0^*  w slots  on  the  average  before 


forwarding  a new  label  to  R^. 


th 


then 


Observe  that  if  tj__^  is  the  expected  end  of  the  ( j — 1 ) cycle 


t . 
J 


+ 0 . w 
3-1 


(2) 
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Let.  q and  q(.  lie  the  probability  of  .successful  transmission  on 
the  repeater  to  station  hop  and  on  the  station  to  repeater  hop,  re- 
spectively , 


OBJECTIVES : 

1.  Determine  q^,  as  a function  of  m,  I,  ft,  w. 

2.  Determine  1,  given  (1). 

ANALYSIS: 

The  phenomenon  is  complicated.  The  How  chart  of  Figure  2 
shows  the  steps  involved  in  initialization. 

Consider  (for  simplicity)  a single  P^  . It  is  broadcasting 
ROP's.  When  a ROP  reaches  S,  the  station  tries  to  send  a label  to 
. This  may  require  several  retransmissions  by  S (in  the  mean- 
time R^  has  issued  more  POP's).  When  R^  finally  receives  a labex 
it  stops  issuing  the  ROP's  and  it  issues  one  ETE  ACK.  If  the  ETE 
ACK  is  received  by  S,  R^  is  eliminated  from  the  queue;  otherwise, 

S retransmits  the  label  until  it  is  received  and  acknowledged  by 

♦ 

R.  . 

l 

2.3  Determination  of  q.  and  q.. 

R 


Let  SD,  Gd  be  the  total  successful  and  total  transmission  rate, 
respectively,  from  the  repeaters  to  the  station.  Let  Sg,  be  the 
total  successful  and  total  transmission  rate,  respectively,  from 
the  station  to  the  repeaters. 

Let  SD  , Gn  be  the  respective  rates  from  R.  to  S.  Let  Sc  D , 

R . R . 1 • 

11  1 

Gc  be  the  respective  rates  from  S to  R. . 
b , R . 1 


2.5 
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We  make  the  following  simplifying  assumptions: 


SR.  SR/m;CR  ■ GR/m;Ss^R  Ssy/m ;(,S , U . “ C\c/In  ( 


3) 


Consider  a specific  slot  in  the  ALOHA  broadcasting  procedure, 
Repeater  R^  is  idle  if  the  following  conditions  are  all  sat- 
isfied. 


1.  The  station  is  not  broadcasting. 

2.  The  neighbors  interfering  with  R,.  are  not  broad- 
casting . 


3.  lb  itself  is  not  broadcasting. 
Thus , 


Prob  (lb  is  idle) 
where , 


(1-0  ) !1 

W 

J l 


( 1_"gr  . ) (1-gr.} 
1 i 


(4) 


V = (R.|d(R.,R.)  < p , \i  = range  of  repeaters } -R . 

3 1 1 1 

As  per  assumption  3 above,  |VR  | --  I,\b;  also  we  postulated  that 

Gr  = GR/m ; thus  1 

1 i 

G 1 + 1 

Prob (R . is  idle)  = (1-G„) (1 — -)  (5) 

jl  o m 

The  station  is  idle  if  1.  and  2.  below  are  satisfied: 
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1 . 

The  s fni. ion 

is  not  broadcasting. 

2. 

No  repeater 

is  broadcasting. 

Thus , 

m G 

Prob  (S  is  idle)  = (1-CQ)  II  (1-Cn  ) - ( 1-G„ ) ( 1- m (6) 

S i=1  R.  S m 

Using  the  above,  we  can  compute  the  expected  throughput  on  the  two 
links. 


Sn  - l (traffic  put  out  by  R.)  Prob  (S  is  not  disturbed 
R i i 

by  any  1 epeater  except  R^  , and  is  idle  itself) 


Gt,  Gt  , Gn  , 

y R n R.m-1 , , , R . m- 1 , . . 

^ m in  S R m S 


(7) 


Ss  = 1'  (traffic  put  out  by  S to  R^)  Prob  (R^  is  not 
disturbed  by  any  device  except  S) 


G G 1 + 1 G 1 + 1 

l = G (1 — ~) 

^ m m S m 


Now 


Rc  = 


S G, 


G T+l 

(1-|) 


K 


qR  G 


R 


(1 — m 1 ( i-G  ) 
lu  b 


(B) 


(9) 
'1(  ) 


If  = p(m)m  and  G^a  monotonically  from  below  as 
Gg  = o (n)m  and  Gg  vb  monotonically  from  below  as 
I = Am,  then 
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(1-p 

(m))m  1 (1-n  (m)m) 

(ID 

(1-p 

. . . A ITH  1 

(m)  ) 

(12) 

The  following  relations  hoi  : 

a.  q_  decreases  as  m increases. 

R 

b.  Lim  qR  = e a ( 1— b ) 
m 

c.  q^  decreases  as  m increases. 

d.  Lim  qs  = ( (e  a) 

(on  the  other  hand,  if  I is  fixed  then  the  limit  is  1) 


Due  to  the  high  complexity  of  the  phenomenon,  we  are  in  a 
position  tc  compute  only  worst  case  analysis. 

A '•O 

Let  G,,  = worst  G„,  GD  = worst  Gn.  The  worst  Gc  is  obtained  when 

o o K K b 

the  station  label  queue  is  assumed  non-empty,  lesultiag: 

6S  = i (13} 

To  obtain  Gn  it  is  assumed  that  all  repeaters  are  transmit  Ling  LOP1^' 
and  that  all  label  packets  to  repeaters  are  successful,  resulting 
in  the  additional  ETE  ACK  from  repeaters.  This  results: 


R 


1 

w 


(14) 
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Substituting  those  values  in  the  expressions  for  C-  and  , we 
derive : 


,,  1 l.m-1.  1. 

f*K  “ (1"n'mw>  (]-w5 


(i-i  ~1) 1+1 

' Si  mw' 


Observe  that  tne  consequence  of 
made  above  is: 


(15) 

(16) 

the  uniform  traffic  assumption 


qR  = qR*  lor  all  i 
i 

/\  A* 

E(Ai)  i E (A) , for  all  i 
E(B^)  E(B),  for  all  i 

Naturally,  not  every  repeater  enters  (leaves)  the  queue  at  the 
same  time,  but  if  the  experiment  were  carried  out  a large  number  of 
times,  on  the  average,  would  enter  (leave)  the  queue  at  E(A)  (E(B)). 
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3. 


1 N LTI ALIZATTON  COMPUTATIONS 


*r- 


Consider  ROP  transmissions  by  R.  with  probability  of  success 
Then  A^  is  a random  variable  geometrically  distributed.  Thus, 


1/q  represents  the  expected  number  of  transmissions  by  R.  before 
R 1 

the  ROP  reaches  S 
(on  the  aver  ::e) , 


the  ROP  reaches  S.  Since  the  ROP  is  put  out  by  R^  every  P slots 


E(A.) 


(17) 


*R 


Consider  label  transmissions  by  S.  1/q^  i epresents  the  ex- 
pected number  of  label  transmissions  by  S.  I f we  knew  the  cycle 
length  at  each  retransmission,  we  could  add  these  values  to  obtdixn 
the  time  necessary  to  deliver  the  label  to  R.. , such  cycle  length 
is  unavailable,  hence,  we  r ist  assume  the  worst  case  of  mw.  Thus, 

— is  the  time  required  to  deliver  one  label  to  R..  However,  it 
qc  i 

1 

takes  - such  receptions  and  ETE  ACK  transmissions  before  the  station 
qR 

receives  an  ACK  and  eliminates  this  repeater  from  the  queue. 

Finally; 


E(B.)  = ^ 

1 qR  q-S  qR 


= i 


(IP) 
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4.  OPTIMIZATION  OF  INITIALIZATION  TIME. 


It-  .is  clear  that,  if  ft  is  small,  the  delay  between  ROP  trans- 
missions is  small,  and  this  would  apparently  speed  up  initializa- 
tion. However,  if  ft  is  small,  the  traffic  genera  tea  by  the  re- 
peaters is  large,  causing  and  to  decrease  and  thus  slow  down 
initiali zation . 

On  the  other  hand,  if  ft  is  large,  the  traffic  is  small  causing 
qD  and  qc  to  increase  and  thus  speeding  initialization;  however, 
since  the  wait  between  ROP  transmissions  is  large,  the  initializa- 
tion time  may  still  increase. 

The  same  holds  true  for  w.  It  is  apparent  that  optimal  values 
of  w and  ft  must  exist.  From  the  previous  section; 


i (ft  , w)  = - 


1 ■ <^+^)J,T1 

ft  mw  w 


ft  + 


11 


mw 

(X-i) 

ft  mw 


It  1 


-- 1 


19) 


Close  form  optimization  of  this  function  of  two  variables  is  un- 
attainable. We  must  thus  resort  to  numerical  techniques,  in  par- 
ticular point-to-point  search,  which  works  well  in  view  of  the 
limited  domain  of  w and  ft  (w  and  ft  must  be  integers  since  we  are 
in  a slotted  environment) . 

The  graphs  of  Figure  3 through  Figure  8 show  the  trends. 

Figure  3 depicts  the  behavior  of  1 as  a function  of  ft  and  w. 
Basically,  we  see  that  the  initialization  time  to  the  right  of  the 
o timum  increases  rather  slowly,  while  to  the  left  of  the  optimum 
increases  drastically;  hence,  we  should  over  estimate  ft*  rather 
than  under  estimate.  Observe  also  that  the  station  rate  is  much 
larger  than  the  repeater  rate;  this  is  so  since  S must  serve  m 
repeaters . 

Figure  4 shows  the  relation  between  ft*  and  w at  I*.  It  demon- 
strates that  if  the  station  transmits  at  a higher  interval,  so  should 
the  repeaters. 
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Figure  5 describes  the  behavior  of  f*  as  a function  of  w, 
pointing  at  the  necessity  of  selecting  the  optimum  value  of  w. 

Figure  6 depicts  the  behavior  of  -F  as  a function  of  m.  'Hie 
graphs  show  that  as  in  increases,  the  optimal  initialization  time 
increases  approximately  linearly  with  m (for  fixed  interference 
level  I ) . 

Figure  7 and  Figure  8 illustrate  the  behavior  of  i as  a func- 
tion of  I.  The  graphs  show  that  as  I increases,  the  optimal  ini- 
tialization time  increases  approximately  linearly  with  I (for  fixed 

m)  . 
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4. 


FIGURE  3:  INITIALIZATION  TIME  AS  FUNCTION  OF  Cl  AND 
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4.4 


FIGURE  4:  RELATION  OF  w AND  Q* , AT  THE  OPTIMUM  INITIALIZATION  TIME 
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4.5 


FIGURE  5:  RELATION  BETWEEN  w AND  ± AT  THE  OPTIMUM  REPEATER  RAT 
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IGURE  7:  RELATION  BETWEEN  I*,  SI*,  w*  AND 
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5 . APPROXIMATE  OPTIMAL  VALUES 

In  lieu  of  a close  form  solution  for  the  optimal  values,  wc 
can  carry  out  a satisfactory  empirical  evaluation  of  the  optima. 
Part  A of  Table  1 gives  the  optima  for  various  values  of  m 

A 

and  I.  Strikingly,  w*  is  constant  over  a large  range  of  m ana  I. 
In  fact  w*  = 4 in  al  1 cases,  except  when  I is  very  close  to  m, 

At 

in  which  case  w*  = 5.  We  thus  conclude  that: 

/A 

w*  = optimal  station  rate,  is  nearly  independent  of  m 

and  I (except  in  case  of  extreme  interference)  and 
w*  = 4 slots. 


Part  B of  the  table  s the  values  for  the  quantity 


, , lxm+l/2+l 

h = ml  + -) 
m 


(20) 


We  see  that  V*  is  very  close  to  this  expression,  we  thus  declare 
that 


= m(l+i)ra+I/2  + 1 
m 


(21) 


Observe  1=2, 


0*  ~ me  (22) 

where  e = 2.718  ...  (more  generally,  this  holds  also  for  the 
case  where  I is  small  and  independent  of  m) . 

A 

Actually,  both  approximations  for  tend  to  overstate  the  value 
of  the  optimum,  as  a comparison  of  Column  B in  Table  1,  with  column 
(d)  reveals.  We  conclude  that: 
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A 

H 

C 

D 

(a) 

(b) 

(c) 

(cl) 

(e) 

Comp  lit  ed 

Computed 

m 

I 

w* 

f * 

Approxi- 

( rounded) 

( rounded ) 

mat ion 

3 

2 

4 

12 

63 

12 

75 

5 

2 

4 

17 

95 

18 

106 

7 

n 

€~ 

4 

22 

128 

23 

138 

9 

2 

4 

27 

161 

28 

1 /I 

(2) 

11 

2 

4 

33 

195 

34 

204 

50 

2 

4 

133 

84  2 

140 

858 

100 

2 

4 

260 

1672 

275 

1698 

1 000 

2 

4 

2 566 

10610 

2722 

16818 

10000 

O 

4 

24626 

165944 

27186 

168018 

100 

10 

4 

273 

1724 

287 

1750 

100 

50 

4 

337 

1986 

350 

2049 

(1) 

100 

90 

5 

429 

2250 

427 

f- 

2431 

11 

4 

4 

36 

208 

37 

219 

11 

6 

4 

39 

221 

40 

237 

(1) 

11 

8 

4 

42 

235 

44 

256 

:•  i 

10 

5 

' 49 

247 

48 

278 

TABLE  1:  NUMERICAL  COMPARISON  OF  EXACT  SOLUTION  TO  APPROXIMATE 

MODEL  AND  APPROXIMATE  SOLUTION  TO  APPROXIMATE  MODEL 
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ft*  = optimal  repeater  rale,  is  an  increasing  function  of 
m and  I,  and  for  small  values  of  I,  u*  me. 

Figure  9 compares  the  approximate  optimal  values  of  the  approxi- 
mate solution  with  the  exact  optima]  values  of  the  approximate 
solution  using  the  results  of  Table  1.  Observe  the  close  fit  for 
a large  range  of  m. 

With  the  above  approximations,  we  can  compute  the  worst  case, 
approximate  initialization  time,  as  follows: 

= (1  - (K^))m'1(1"i)  (23> 


(_Lf  Jo  im_1  ( 1-— ) 

'me  4m' 1 ' 4' 


3 f , 4-le , m- 1 _ 2 r i J 4-10)740,  m-  1 

4 1 4em  4 1 1 m J 


2 (4+e)/4e 

4 e 


Similarly , 


(_i+  i) ]i+l  = [i 

me  4m 


(4+e) / 4e, m 

j m 

m 


(4+e)/4e^  (I+l)/m  ^ ^ 54)(I+l)/m 


Finally, 


~ + ft  , mw 

■%  m 


* 40  .4  (.  54)(lTTT7m' 


= 6 . 79m  + 


(.54) 


10m 

(I+l)/m 
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MODEL  AND  APPROXIMATE  SOLUTION  TO  APPROXIMATE  MODE 
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I art  C oi  Table  1 shows  for  comparison  the  computed  values  of  T*. 

So  far  i has  b^on  quoted  in  a number  of  slots  in  the  ALOHA 
procedure.  Typical  value  for  slot  length  for  an  ROP  is  5xl0~“* 
seconds.  Hence,  for  example,  in  a 100  repeater  system  at  inter- 
ference level  2,  it  would  take  0.84  seconds  to  initialize  a typical 
repeater . 
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6 , CONCLUSION 

In  this  chapter  a worst  case  study  for  initialization  has  been 
conducted;  it  is  believed  that  the?  same  trends  with  respect  to  the 
dependency  of  f on  I and  n,t>  etc.,  hold  true  in  the  actual  initial- 
ization procedure.  To  ascertain  this  an  exact  model  via  a Markov 
chain  is  being  developed,  and  will  be  compared  to  the  model  described 
presently . 

Other  approximate  models  (multiple  retransmissions,  finite 
buffer  size,  2-hops)  are  under  investigation. 
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